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ON  THE  POWER  FUNCTIONS  OF  TEST  STATISTICS 


IN  ORDER  RESTRICTED  INFERENCE 


Hari  Mukerjee,  Tim  Robertson,  and  F.  T.  Wright 


SUMMARY 


"We  study  the  power  functions  of  both  the  likelihood  ratio  and  con¬ 
trast  statistics  for  detecting  a  totally  ordered  trend  in  a  collection  of 
means  of  normal  populations.  Monotonicity  properties  are  found  and  both 
radial  limits  and  limits  along  lines  parallel  to  the  cone  of  points  satis¬ 
fying  the  trend  are  examined.  An  optimal  contrast  test  for  testing  a  trend 
as  a  null  hypothesis  is  derived. 


AMS  1980  subject  classifications:  Primary  62F03;  Secondary  62F04. 


Key  words  and  phrases:  Order  restricted  tests,  isotonic  inference,  power, 
likelihood  ratio  tests,  contrast  tests. 


INTRODUCTION.  We  consider  the  powers  of  statistical  tests  for  detecting 
a  trend  in  a  collection  of  population  parameters.  The  statistical  litera¬ 
ture  contains  a  number  of  such  tests  and  a  detailed  summary  of  this  research 
up  to  about  1971  is  given  in  Barlow  et  al .  (1972).  More  recent  summaries 
are  given  in  Bartholomew  (1983),  Lee  (1983)  and  Robertson  (1984).  We 

restrict  our  attention  to  the  case  in  which  the  parameters,  ,u2 . 

are  the  means  of  normal  populations  and  the  trend  restriction  requires  them 
to  be  totally  ordered.  To  be  specific  we  consider  the  trend  H-j :  s  u2 

£  • • •  £  |i^.  One  approach  to  detecting  such  a  trend  is  to  test  homogeneity, 
HQ:  =  u2  =  •••  =  nk»  versus  i.e.,  H]  holds  with  ^  <  uk 

(cf.  Bartholomew  (1959  a,b;  1961)).  Even  in  the  case  of  normal  means  the 
results  concerning  the  powers  of  these  restricted  tests  are  very  exiguous 
and  consist  primarily  of  comparisons  with  the  powers  of  other  tests,  such 
as  the  unrestricted  tests  of  HQ  versus  t  for  some  i  /  j. 

In  fact,  as  far  as  we  can  determine,  the  first  mention  of  the  fact  that 
Bartholomew's  tests  are  unbiased  occurs  in  Robertson  and  Wright  (1982). 

The  biases  of  other  restricted  tests  are  examined  in  Dykstra  and  Robertson 
(1983). 

Assuming  independent  random  samples  from  normal  populations, 
Bartholomew  (1959  a,b;  1961)  studied  the  likelihood  ratio  tests  (LRTs) 
for  Hq  versus  H-j-Hq  assuming  in  one  case  that  the  population  variances 
are  known  and  in  the  other  that  they  are  unknown  but  equal  (partially 
ordered  trends  were  also  considered).  We  focus  attention  on  the  case  of 
known  variances.  Results  concerning  the  unknown  variances  case  follow  by 
conditioning  arguments  in  the  last  section.  Implementation  of  Bartholomew's 
test  procedures  can  be  difficult  for  k  >  5  if  the  so-called  weights  are 


not  equal.  (For  the  known  variances  case,  the  weights  are  the  precisions, 

2 

n  of  the  sample  means  as  estimates  of  the  population  means.)  This 

difficulty  is  mainly  due  to  the  fact  that  the  level  probabilities  involved 
in  the  null  hypothesis  distribution  of  the  test  statistic  are  extremely 
difficult  to  compute  in  such  cases.  This  theory  is  discussed  at  length  in 
Chapter  3  of  Barlow  et  al.  (1972).  Robertson  and  Wright  (1983)  proposed  an 
approximation  for  the  level  probabilities  for  the  case  of  total  order  and 
Pillers  et  al.  (1984)  gives  a  computer  routine  for  implementing  this  approx 
imation. 

Partly  because  of  the  difficulty  involved  in  applying  Bartholomew's 
procedures,  several  researchers,  including  Abelson  and  Tukey  (1963)  and 
Schaafsma  and  Smid  (1966),  considered  tests  based  upon  contrasts  (cf.  Sec¬ 
tion  4.2  of  Barlow  et  al .  (1972)).  Denoting  the  sample  means  by  X\ , 

1  <:  i  s  k,  these  contrast  tests  are  based  upon  statistics  of  the  form 
|<  _ 

Tc  =  Ei=1  c^.  where  c  =  (c^ ,c2»* • • ,ck)  is  a  vector  of  predetermined 
constants  (E.^  ci  =0).  One  attraction  of  these  constrasts  is  the  fact 
that  their  distribution  at  any  point  n  =  (n-j  ,u2» •  •  •  ,uk)  €  R  is  normal. 
With  Hq,v  6  R  fixed  (E^_-|  =0),  the  uniformly  most  powerful  (UMP) 

test  of  Hq:  \jl  =  iuq  versus  :  m  €  [|Uq  +  bv;  b  >  0]  rejects  Hq  for 
large  values  of  T^  (use  the  Neyman-Pearson  Theorem  and  note  that  this 
test  is  UMP  for  fixed  b  >  0).  Since  this  test  does  not  depend  on  Hq€Hq  i 
is  not  surprising  that  contrast  tests  are  very  powerful  in  some  subregion 
of  the  alternative.  However,  even  for  moderate  k,  there  are  other  sub- 
regions  of  the  alternative  where  the  power  of  the  contrast  test  does  not 
compare  favorably  with  the  power  of  Bartholomew's  LRT.  While  the  LRT  is 
not  most  powerful  at  any  particular  point,  it  does  maintain  a  more 
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uniformly  reasonable  power  over  all  of  H-j .  One  explanation  of  this  fact 
is  given  in  Theorem  4.3  of  Barlow  et  al.  (1972),  which  can  be  interpreted 
to  say  that  the  LRT  is  based  on  an  "adaptive"  contrast  statistic.  In  other 
words,  the  parameters  are  estimated  from  the  data  and  then  the  contrast 
coefficients  are  chosen  so  that  the  test  has  a  relatively  high  power  at  the 
estimated  point. 

In  the  references  cited  earlier,  Abelson  and  Tukey  (1963)  and 
Schaafsma  and  Smid  (1966)  derived  optimal  contrast  tests.  In  the  former 
they  obtained  the  contrast  coefficients  which  maximize  the  minimum  power 
over  all  points  equidistant  from  Hq  and  in  the  latter  those  that  minimize 
the  maximum  loss  of  power  as  compared  to  the  most  powerful  test  in  a 
restricted  class  of  procedures.  The  powers  of  these  optimal  contrast  tests 
are  compared  with  that  of  the  LRT  in  Section  4.2  of  Barlow  et  al .  (1972). 
Their  conclusion  is  that,  while  it  is  difficult  to  improve  on  the  LRT,  for 
a  total  order  and  small  k  these  contrast  tests  provide  viable  alterna¬ 
tives  to  the  LRT. 

In  the  above  approach  to  detecting  a  trend  the  hypothesis  of  homo¬ 
geneity  is  a  "dummy"  hypothesis.  Control  of  the  a  level  provides  pro¬ 
tection  against  confirming  the  trend  when  it  is  not  present.  Robertson  and 
Wegman  (1978)  and  Robertson  (1978)  studied  LR  procedures  for  testing  a 
trend  as  a  null  hypothesis.  The  a-level  of  these  tests  controls  the 
probability  of  denying  the  trend  when  it  is  present.  The  null  hypothesis 
distributions  of  the  LRT  statistics  involve  the  same  level  probabilities 
that  complicate  the  use  of  Bartholomew's  test.  In  Section  4  we  derive  an 
optimal  contrast  test  for  H.  versus  H0:  u.  >  n- ,,  for  some  i. 
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In  Section  3,  we  study  the  power  functions  of  the  LRTs  for  testing  Hq 
versus  H-j-Hq  and  for  versus  ^  The  power  functions  for  the  latter 
tests  are  more  complicated  and,  in  a  sense,  more  interesting.  Some  prelim¬ 
inary  results  for  studying  the  powers  of  the  LRTs  are  developed  in  Section 
2.  The  competing  contrast  tests  are  discussed  in  Section  4. 

Throughout  this  paper  "X ^  » *  *  * »\  denote  the  sample  means  of  inde- 

r\  XL 

pendent  random  samples  with  X.  ~  n^.  ,o^/n^ ) ,  n^  the  size  of  the  i— 

2  f  k 

sample,  and  a f  the  variance  of  the  i—  population.  Assume  that  the 

2 

variances  are  known  and  set  w.  =  n ./a.  for  i  =  l,2,***,k.  Let  H„  and 

iri  0 

H-j  denote  the  hypotheses  specified  earlier  as  well  as  the  corresponding  sub- 
1/ 

sets  in  R  .  The  set,  Hq,  is  a  subspace  and  is  a  closed,  convex 

1/ 

cone.  Let  (•,*)u  denote  the  inner  product  on  R  defined  by 
(x,y)w  =  w.x.y.  and  let  ||*||w  denote  the  corresponding  norm. 
Bartholomew's  test  of  HQ  versus  H-|-Hq  rejects  HQ  for  large  values  of 

(1.1)  T0]  =  -2  In  A  ■  wf(5ru)2  =  IliT-ueJ!2 

a  k  — •  k 

where  A  denotes  the  likelihood  ratio,  u  =  E.  ,  w.X./E.  .  w. ,  e.  is  a 

1 =  I  1 =  I  1  K 

k-dimensional  vector  of  ones  and  j I  =  (u^ » • • * )  minimizes 

£j_-|  w^(g..-X..)  subject  to  g  €  H^.  In  other  words,  u  is  the  projection 

of  J  =  (X^ , Xg , * • • ,Xk)  onto  with  respect  to  the  distance  d(x,y) 

=  j|x-yi|w.  We  will  also  denote  this  projection  by  EW(X"I  H-j ) .  With  the 

y 

closed,  convex  cone  D  {x:Z^_i  x.w,.  =  0]  denoted  by  Cq^  ,  Theorem  1.5 
of  Barlow  et  al.  (1972)  can  be  used  to  show  that  E^(T|  Cq-j  )  =  EWCX|  )  -  u 

_  O 

and  so  TQ^  =  j|Ew(X  (  C01 ) ^W*  Therefore>  an  acceptance  region  for  TQ1 
be  written  as  [x  €  :  ||Ew(x  [  Cq-j ) j]^  s;  t}  for  some  t  >  0. 

The  LRT  of  versus  H ^  rejects  for  large  values  of 


can 
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0-2)  T)2  =  -2  In  A  =  2*.,  w^X,-^)2"  ||X-EW(X  |  H,)U2. 

If  C-|2  denotes  the  dual  of  H-j  (also  called  the  polar  or  conjugate  of 

H.| ) ,  which  is  a  closed  convex  cone  whose  definition  is  given  in  the  next 

section,  then  Theorem  1.5  of  Barlow  et  al.  (1972)  shows  that  E^(x  |  C^) 

=  x-Ew(xi  H-j ) ,  and  so  =  !|EW(X  i  The  acceptance  regions  for 

T-j 2  are  of  the  form  {x  €  :  [jE^(x  j  C^)!]^  *  t]  with  t  >  0. 

In  Section  3  we  consider  the  question  of  unbiasedness  for  Tg^  and 

T-j 2  as  well  as  the  radial  behavior  of  their  power  functions,  that  is, 

their  behavior  on  the  sets  { 6j_l ;  6^0}  for  various  n.  The  behavior  of 

these  power  functions  in  other  directions  is  also  discussed.  Robertson  and 

1/ 

Wright  (1982)  considered  the  relation  on  R  ,  x<y  provided  y-x  6  . 

They  showed  that  T  ^  and  its  power  function  are  isotone,  and  T12  and  its 
power  function  are  anti  tone  with  respect  to  <.  This  implies  that  if 
u  €  H-j ,  then  the  power  of  TQ1  (T12)  is  nondecreasing  (nonincreasing)  in 
6  on  :  -<*  <  6  <  *}.  Their  results  concerning  the  stronger  relation 
«  show  that  the  power  of  T^  is  monotone  in  a  larger  set  of  directions, 
but  these  techniques  do  not  apply  to  T^2  because  it  is  not  antitone  with 
respect  to  «.  However,  because  of  the  strong  similarities  in  the  accept¬ 
ance  regions  for  the  tests  Tg^  and  T^,  one  might  conjecture  that  T12 
is  also  monotone  in  this  larger  set  of  directions.  Using  the  geometric 
arguments  of  Section  2,  this  is  shown  to  be  the  case. 


2.  SOME  PROBABILITY  INEQUALITIES.  The  probability  inequalities  derived 
in  this  section  will  be  used  in  the  discussion  of  the  montonicity  of  the 
power  functions  of  Tq-|  and  T^.  Using  the  techniques  in  Bartholomew 
(1961),  one  can,  at  least  in  principle,  obtain  analytic  expressions  for 
these  power  functions  in  terms  of  several  multiple  integrals,  but,  as  func¬ 
tions  of  the  distance  of  a  mean  vector  from  Hg  and  its  direction,  they 
are  extremely  intractable  even  in  the  case  of  equal  weights  and  k  =  3. 

(See  Section  3  for  further  discussion.)  We  have  resorted  to  geometric 
arguments  involving  integrals  of  symmetric  unimodal  densities  over  convex 
sets  that  have  certain  symmetry  properties.  In  a  sense,  the  results 
obtained  are  generalizations  of  Anderson's  (1955)  work  on  similar  integrals 
over  symmetric  (about  the  origin)  convex  sets,  but  in  the  present  work,  the 
statistics  are  projections  on  closed,  convex  cones  which  are  not  subspaces, 
and  thus  the  sets  involved  have  only  partial  symmetry. 

The  basic  idea  is  the  following.  Let  P  (•)  denote  the  n(n,I) 

k  k 

probability  distribution  on  R  for  each  n  €  R  ,  and  note  that 
P^(A)  =  Pg(A-u)  for  all  measurable  AcR  ,  If  A  is  the  acceptance 
region  for  one  of  the  tests  considered  here,  then  A  is  closed  and  convex. 
If  S  is  the  subspace  generated  by  a  mean  vector  n  and  S'1  is  the 

k 

orthogonal  complement  of  in  R  ,  then,  for  some  directions  of  u, 

it  is  possible  to  decompose  A  into  disjoint  subsets  A'  and  A",  with 

A'  a  closed,  convex  set  symmetric  about  Sx  and  A"  on  one  side  of  Sx. 

iu  u 

This  will  enable  us  to  prove  the  monotonicity  of  P^(A)  =  Pq(A-6h)  in 
6  is  0  for  such  directions. 

Because  we  anticipate  the  application  of  the  results  of  this  section 
to  more  general  types  of  cones  than  those  considered  in  this  paper,  and 


because  we  believe  the  results  concerning  projections  on  closed,  convex 
cones  in  a  real  Hilbert  space  are  of  interest  in  themselves,  we  consider  a 
more  general  framework  than  is  needed  for  this  paper. 

Let  H  denote  a  real  Hilbert  space  with  inner  product  (•,•)  and 
norm  [[ - 1| .  If  C  is  a  closed,  convex  cone  in  H  and  x  €  H,  then  E(x  |  C) 
will  denote  the  unique  projection  of  x  onto  C,  i.e.,  E(x|  C)  is  the 
unique  element  of  C  which  minimizes  |[x-y[|  as  y  ranges  over  C.  Theorem 
7.8  of  Barlow  et  al.  (1972)  characterizes  E (x  J  C)  as  follows: 

(2.1)  E(x[C)€C,  (x-E(x  |  C) ,  E(x  [  C) )  =  0,  and  (x-E(x|  C),y)  s  0 
for  all  y  €  C. 

It  follows  from  (2.1)  that  E(ax  (  C)  =  aE(x[  C)  for  a  ^  0  and  x  €  H, 
and  that 

(2.2)  (x,E(xi  O)  =  (E(x  I  C),E(x|  C))  =  |jE(x  [  C)||2  for  x  6  H. 

The  dual  of  C,  which  is  denoted  by  C*,  is  defined  by 

C*  =  (x  €  H  :  (x,y)  s:  0  for  all  y  6  C}.  Clearly,  C*  is  a  closed,  convex 

cone  and  using  the  definition  of  C*,  (2.1),  and  (2.2),  we  see  that 

(2.3)  E ( x  J  C* )  =  x  -  E ( x  (  C )  and  ||E(x  [  C*)[|2  =  [[x|[2  -  |[E(x  (  C)[|2. 

In  the  Appendix  it  is  shown  that 

(2.4)  (C*)*  =  C. 

Throughout  this  section  C  will  denote  a  closed,  convex  cone  in  H, 

C*  its  dual,  and  A  the  closed  set  {x  €  H:  |[E(x  |  C)(|  £  t}  for  some  t  >  0 
For  ij  €  H  let  =  £bp:  b  £  0),  S  =  {b^i:  -®  <  b  <  ®] ,  and  let  be 


the  orthogonal  complement  of  S  . 

We  repeat  for  the  reader's  convenience  a  result  from  Robertson  and 
Wegman  (1978)  which  will  be  used  several  times  in  the  sequel. 

Lemma  2.1.  If  x  €  C,  then  for  any  y  €  H,  [|x+y-  E(x+y  j  C)[| 

^  !|y- E(y  I  C)||  or,  equivalently,  ||E(x+y  I  C*)|j  *  ||E(y  |  C*)||. 

Next,  we  state  several  lemmas  which  are  proved  in  the  Appendix.  Let 
-D  =  {-x  :  x  €  D}. 

Lemma  2.2.  Suppose  pt  €  C  U  (-C*).  For  x  €  H  and  0  s  b  s  2, 

(2.5)  I|E(x- bE(x-px0  J  C^)  J  C)i|  s  ltE(x|C)U  for  all  pQ6Sj. 

So,  if  (pi,piQ)  =  0  and  y  €  A-piQ,  then  y-bE(y(  C^)  €  A-piQ  for  0sbs2. 
Lemma  2.3.  For  x,y  €  H, 

!iE(x+y(  C)||  *  || E (x  |  C)+E(y  |  C)U  s  ||E(x  [  C)|J+[|E(yl  C)||. 

Using  Lemma  2.3,  we  see  that  A  is  convex. 

Lemma  2.4.  Let  S  be  a  closed  subspace  of  H. 

(a)  If  Sc:C,  then  (x-E(x  }  C),v)  =  0  and  E(x-v  [  C)  =  E(x  [  C)  -  v  for 

all  x  €  H  and  v  €  S. 

(b)  If  See,  then  E(E(x  |  C)  [  S)  =  E(x  |  S)  for  all  x  €  H. 

(c)  If  See,  then  E(x  j  C)  -  E(x  {  S)  =  E(x  [  C  fi  Sx)  for  all  x  €  H. 

(d)  If  CcS,  then  E(E(x  [  S)  |  C)  =  E(x  |  C)  for  all  x  €  H. 

The  next  result  identifies  the  portion  of  A  that  is  symmetric  about 
SO  Define  A+  to  be  {x  €  A :  E(x  (  S^)  =  by  for  some  b  s  0]  and  B  to 


be  [x -  bE(x  j  S^)  :  x  €  A  ,  Os  bs  2}. 

Theorem  2.5.  If  p  €  C  U  (-C*)  and  is  any  vector  in  H  with 

(u,un)  =  0,  then  B-nn  c  A-nn,  B-|jn  is  symmetric  about  SA,  i.e.. 


x  €  B-un  implies  x-  2E(xj  S  )  €  B-un,  and  B-nn  is  closed  and 


convex. 


Proof.  First  assume  the  theorem  is  true  with  =  0.  Then  the  first 
and  the  last  conclusions  follow  immediately  for  arbitrary  Hq  6  H.  The 
second  conclusion  follows  when  =  0  by  noting  that,  when  B  is 

symmetric  about  SA,  x  €  B-^q  implies  x+vig-  2E(x+uq[  S  ) 

=  x+Hg-  2E(x  |  S^)  €  B,  so  that  x-  2E(x  j  S  )  6  B-Uq.  Thus  it  is  suffir  t 
to  prove  the  theorem  for  Hq  =  0. 

For  x  €  A+,  E(x  |  C  )  =  E(x  |  S  )  and  so  BcA  by  Lemma  2.2.  If 
x  €  B,  then  there  exists  y  €  A+  and  b  6  [0,2]  with  x  =  y-  bE(y  |  Sy). 
Thus  by  the  linearity  of  E(-  J  S  ), 

x-2E(xiSu)  =  y-bE(y|S(j)-2(l-b)E(y|Sia)  =  y-  (2-b)E(y|S^  6  B. 

Next,  we  show  that  B  is  closed.  Since  S  is  1-dimensional  it  is  closed, 

u 

and  so  is  SA.  Because  A+  is  the  intersection  of  the  closed,  convex  sets 

u> 

A  and  {x  €  H:  E(x  |  S  )  =  bja,  b  s  0],  it  is  closed  and  convex  also.  Now 
A+  =  [x -  bE(x  J  S^) :  x  €  A+,  0  s  b  s  1],  and  if  we  define 

A”  =  {x-  bE(x  [  S  ):  x  €  A+,  1  s  b  £  2],  then  A"  is  the  reflection  of  A+ 

across  SA.  It  is  then  easily  verified  that  A"  is  closed  and  convex. 
Hence,  B  =  A+  U  A"  is  closed. 

Since  A+  and  A"  are  both  convex,  to  show  that  B  is  convex,  we 

need  only  show  that  ax+  (l-a)y  €  B  for  0  <  a  <  1,  x6A+H  (A")c  and 

y  €  A"  fl  (A+)c.  Now  y-  2E(y  J  S .)  €  B  and  E(y- 2E(y|S  .)|S  .)  =  -E(y|S  ) 


implies  that  y-  2E(y  [  S^)  €  A  .  Thus,  v  =  ax  +  ( 1  -a ) (y  -  2E(y  |  S^))  €  A+ 

and  v-  bE(v  [  S  )  €  B  for  0  s  b  s  2.  Now,  E(x  [  S  )  =  c  u  and  E(y  [  S  ) 

=  cya  with  cx  >  0  and  <  0.  Then  E(v[  S^)  =  (acx~  (l-a)cy)pi  and 

v-  b(ac  -  (l-a)c  )u  €  B  for  all  0sbs2.  But  0  <  -2(l-a)c 

x  y  y 

<  2(ac  -  (l-a)c  ) .  Thus 

A  y 

ax+  (l-a)y  =  v+  2  ( 1  -  a )  cyj.  =  v-[(-2(l-a)cy)/(acx-(l-a)cy)]E(vlS^)  €  B 
and  the  proof  is  completed. 

k 

For  the  remainder  of  this  section,  we  take  H  to  be  R  and  at  first 

k 

consider  the  usual  inner  product.  For  pi  6  R  we  let  P^  denote  the 
n(y,I)  probability  distribution  on  Rk. 


k  k 

Theorem  2.6.  If  pi  €  R  and  E  is  a  closed,  convex  set  in  R  which 


is  symmetric  about  S^,  then  P^(E)  is  nonincreasing  in  &  for  6*0. 


Proof.  We  may  assume  that  n  f  0  and  [jptjj  =  1.  Introduce  in  Rk  an 


orthonormal  coordinate  system,  so  that  if  x  €  RK  has  coordinates 


(x-|  jXg,  •  •  •  jX^),  then  E (x  [  S  )  =  (x-j  ,0,  •  •  •  ,0) .  Hence,  for  any  measurable 
k 


DcR' 


P0(D>  =  J(J  J  ni=2  cp(x.)dxi)cp(x1)dx1  =  JgD(x1)cp(x1)dx1 


D(x1 ) 


where  D(x^)  =  [ (x^, •  •  •  ,xk)  :  (x^ .x^, •  •  •  ,x^)  €  D}  is  the  x-j-section  of  D 
and  cp  is  the  standard  normal  density.  From  the  symmetry  assumption  on 
E,  E(-x^)  =  E(x^),  and  because  E  is  convex  and  symmetric,  E(ax^) 
c;  E(a'x^)  for  0  £  a'  £  a.  So,  g^  is  symmetric,  nonnegative,  and  non¬ 
decreasing  on  ( -°°,0] ,  and 


V 


%  ■ 


v 


V  ■  L  ■ 

-  _  ^  « 
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P^E)  =  Pq(E-6h)  =  J‘gE(x1)cp(x1-6)dx1  =  JgE(x1+6)cp(x1  )dxr 
Suppose  0  s  6  s  6' .  Let 

I  =  Pq(E-6u)  -  P0(E-6V)  =  J[gE(x-,+6)  -  gE(x1+6/ )]cp(x1)dxr 

Let  c  =  (6+6' )/2  and  b  =  (6'-6)/2  and  note  that  c  a  b  ^  0.  Hence, 

1  =  jCgE(x-|+c-b)  -  gE(x1+c+b)]cp(x1)dx1 

=  J[gE(y-b)  -  gE(y+b)]cp(y-c)dy  =  /h(y)cp(y-c)dy, 

where  h(y)  =  gE(y-b)  -  gE(y+b)  is  antisymmetric  and  h(y)  £  0  for  y  &  0. 
Thus, 

ro 

I  =  S  h(y)[cp{y-c)  -  cp(y+c)]dy. 

0 

2  2 

For  y  0  and  cs  0,  (y+c)  2s  (y-c)  and  so  1  2  0.  The  proof  is 

completed. 

Remark.  Replacing  u  by  -u  in  Theorem  2.6,  we  see  that  P  c  (E) 

“CHi 

is  nonincreasing  in  6  2*  0  also. 

k  k 

Theorem  2.7.  If  n  €  R  and  D  is  a  measurable  subset  of  R  with 

E(x  [  C  )  =  0  for  each  x  €  D,  then  P^(D)  is  nonincreasing  in  6  2  0. 

Proof.  We  assume  f  0  and  ||u|l  =  1.  We  use  the  coordinate  system 

and  notation  used  in  the  proof  of  Theorem  2.6. 

Because  of  the  hypothesis  on  D,  D(x^)  =  <f>  and  9q(xi)  =  0  for  all 

x-j  >  0 .  If  0  s:  6  £  6' ,  then  PQ(D-6ii)  -  Pq(D-6'h) 

=  /  gn(x1)[cp(x,-6) -<p(x,-6'  )]dxr  Since  g^x^O  and 


[1] 


cp( x-j -6 )  *  cp(x^-6/)  for  x-|  £  0  and  0  s  5  s  a' ,  the  last  integral  is 
nonnegative.  The  proof  is  completed. 

If  D  and  O'  are  subsets  of  a  linear  space,  we  write  DSD'  to 
denote  the  direct  sum,  i.e.,  D  ©  O'  =  {x+y  :  x  €  D,  y  6  D'). 

k  I  k 

Theorem  2.8.  Let  pi  €  R  ,  Uq  €  S  ,  and  let  S  be  a  subspace  in  R 
containing  C.  If  n  €  (C  ©  S1)  U  (-C*),  then  ^(A-Hq)  is  nonincreasing 
in  6*0,  and,  if  n  6  C*.  then  P^(A-pig)  is  nondecreasing  in  6*0. 

Proof.  First,  consider  the  case  u  6  C  U  (-C*).  Combining  Theorems 
2.5  and  2.6,  we  see  that  P^(B-Uq)  is  nonincreasing  in  6*0.  If 
x  €  A  n  Bc-uq  then  E(x  [  C  )  =  0,  and,  applying  Theorem  2.7, 

P6p(An  BC-u0)  is  nonincreasing  in  6*0.  With  n  €  C  U  (-C*)  and 
v  €  S'1,  p<5(u+v)(A“Fio)  =  VA^0_6v) *  The  first  conclusion  will  be  estab¬ 

lished  by  showing  A-6v  =  A.  Now,  A-6v  =  C*-6v  :  x  €  A] 

=  {x-6v:  j|E(x  1  C) |I  t}  =  £y  :  j[E(y+6v|  C)||  £  t}.  Applying  part  (d)  of 
Lemma  2.4,  the  linearity  of  E(-  j  S),  and  the  fact  that  E(6v  [  S)  =  0,  we 
see  that 

E(y+6v  I  C)  =  E(E(y+6v  |  S)  |  C)  =  E(E(y  [  S)  I  C)  =  E(y(  C). 

So  A-&v  =  {y  :  IJE(y  |  C) U  s:  t]  =  A. 

For  the  second  conclusion,  we  assume  that  n  €  C*  and  0  s  6  s  6' . 
Since  (6'-6)n  f.  C*,  by  Lemma  2.1,  |[E(x+6'ia|  C ) [j  £  ||E(x+6u  I  C)JJ.  Thus 

x+  6u  €  A  implies  that  x+  6/u  €  A  or  A+6'ucA-6pi,  so  that 
A-  u0-  6#u  c  A-  Hq-  6m .  Hence,  P^  (A-iaQ)  =  pq(A-Uq-6V)  *  PQ(A-u0-6u) 

=  P.  (A-yn).  The  proof  is  completed. 


Theorem  2.8  will  be  used  to  study  the  monotonicity  of  the  power  func¬ 
tions  of  TQ1  and  in  the  case  of  equal  weights,  i.e.,  w-j  =  w^  = 

•••  =  w^.  The  analogous  results  for  unequal  weights  will  be  established 
next.  Let  w..  >  0  for  i  =  1,2, and  let  W  be  the  kxk  diagonal 

matrix  with  W..  =  w.,  i  =  !,•••, k.  Consider  the  inner  product, 

111  w 

k 

and  norm,  u’liw»  on  R  which  were  defined  in  the  Introduction.  For  C 

a  closed,  convex  cone  in  Rk  (closed  in  the  topology  induced  by  j|*|L), 

let  C  denote  the  dual  of  C  and  Ew(*  |  C)  denote  the  projection  onto 

C  with  respect  to  (•,•).,;  for  S  a  subspace  in  Rk,  let  SxW  denote 

its  orthogonal  complement  with  respect  to  ( - , * ) w ;  for  fixed  t  >  0,  set 

Aw  =  (x  6  Rk  :  JjEw(x  [  C)[[^  s  t];  and  for  u  €  Rk,  let  w  denote  the 

k 

normal  probability  distribution  on  R  with  mean  \jl  and  covariance  matrix 
W"1.  If  W  =  I,  then  we  will  omit  the  subscript  or  superscript,  except 
for  emphasis. 

We  now  establish  two  identities  involving  dual  cones  and  projections 
with  respect  to  (*,Oj  and  (*,*)^*  They  can  easily  be  generalized  to 
the  case  in  which  W  is  a  positive  definite  matrix  or  an  invertible  posi¬ 
tive  operator  on  a  real  Hilbert  space.  In  either  case,  the  inner  product 
(x >y -  (x , Wy ) j . 

1  /2 

Let  W  denote  the  unique  positive  square  root  of  W  and  let  0  be 
any  kxk  orthogonal  matrix.  The  matrix  0  plays  no  essential  role,  how¬ 
ever,  a  judicious  choice,  such  as  a  generalized  Helmert  transformation  in 

the  case  of  a  totally  ordered  trend,  may  help  identify  the  transformed 

1  /2 

parameters  and  visualize  the  transformed  cone.  Let  F  =  OW  and  note 

k  • 

that  F  is  invertible.  For  x,y  €  R  ,  it  is  easily  verified  that 
(x,y)w  *  (Fx,Fy),  |!x|Sw  =  ||Fx||,  and  (x,y)w/(||xj[w|[y|lw)  =  ( Fx , Fy )/ ( 'J Fx|| |[ Fy j; ) , 
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with  Fx  f  0  and  Fy  f  0  in  the  latter  case.  Since  F  is  linear  and 
invertible  FC  is  also  a  closed,  convex  cone  and  C  =  F~^(FC).  The  fol¬ 
lowing  lemma  is  proved  in  the  Appendix. 

Lemma  2.9.  FC*W  =  (FC)*1  and  E(Fx  [  FC)  =  FEw(x  I  C)  for  all  x  €  Rk 

We  now  prove  the  following  generalization  of  Theorem  2.8: 

k  x 

Theorem  2.10.  Let  pi  €  R  ,  €  S  ,  and  let  S  be  a  subspace  in 

Rk  containing  C.  If  u  €  (C  ©  SJ,W)  U  (-C*W),  then  P^  ^AW~W(P  is  non 
increasing  in  6^0,  and,  if  u  €  C*^,  then  vJ^AW-1J0^  1S  nondecreas 
ing  in  6  at  0. 

Proof.  If  X  is  distributed  as  P  u  then  Y  =  FX  has  a  ri(Fu,I) 

U»w 

distribution.  So  P^yW  =  V»,I(F(Vu0»  -  VUjI(  FA,,-^) ,  And 
we  will  apply  Theorem  2.8  to  the  latter  term  with  u,  y  ,  S,  C,  and  A 
replaced  by  Fy,  FnQ,  FS,  FC,  and  FAW,  respectively.  Note  that 
(Fu,Fu0)  =  (u,u0)w  =  0,  FC  c  FS,  (FC  ©  (FS)xI)  U  (-(FC)*1) 

=  F[(C  ©  SaW)  U  (-C*W)]  and  so  pi  €  (C  ©  SxW)  U  (-C*W)  implies 
Fn  €  (FC  ©  (FS)xI)  U  (-(FC)**).  Now,  using  Lemma  2.9, 

FAW  =  tFx  :  |[Ew(x  {  C)||w  s  t]  =  {Fx  :  |[FEw(x  |  C)UX  *  t) 

=  {Fx:  [iE(Fx(  FC)H  *  t} 

Thus,  by  Theorem  2.8,  P^  ^AW~w0'  =  P6F(j  j^A^-F^)  is  nonincreasing 
in  6  it  0. 

The  proof  of  the  generalization  of  the  second  conclusion  of  Theorem 
2.8,  which  is  similar,  is  omitted. 
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3.  LIKELIHOOD  RATIO  TEST.  We  now  apply  the  results  in  Section  2  to  study 
the  power  functions  of  the  LRTs  of  Hq  versus  ^-Hq  and  of  H-j  versus 
H^-  The  power  function  of  T^  is  examined  first.  With  t>0  fixed, 
we  denote  this  power  function  by  tTq^(*)* 

If  k  =  2,  then  the  testing  situation  is  the  classical  one-sided  test 
of  Hi  =  h2  versus  <  u2.  this  case,  rejecting  =  h2  for  large 

values  of  is  equivalent  to  rejecting  for  large  values  of  "X 2  ^  . 

It  is  well  known  that  such  a  test  is  UMP  and  its  power  function  is  increas¬ 
ing  in  h2  -  H]  (cf.  Problem  3.2  on  p.  117,  Lehmann  (1959)).  The  complexity 
of  the  situation  increases  rapidly  in  k.  For  k  =  3  and  Wi  =  w2  =  w^, 
Bartholomew  (1961)  derived  an  expression  for  the  power  function  of  Tqi  . 

(The  derivation  is  also  given  in  Section  3.4  of  Barlow  et  al .  (1972).)  Let 
u  €  R3,  A  =  (I3=1(H.-E(ni  H0))2)172  and  let  3  be  defined  by 


(n2-Ul)A/2  =  A  sin  3  and  (2h3~h2-hi  )/JG  -  A  cos  3. 


The  restriction  n  €  ^  is  equivalent  to  0  s  3  &  tt/3.  With  $  the  stan 
dard  normal  distribution  function  arid  iff(x,t)  =  (x$(xVt)  +  0(xVt) )/0(x) , 


(3.1)  ~  ^Tqi 


>t]  =  6XP("^a2)  j^/2+  3 

tt/6  +  3 


i|i(A  sin  0,t)d6 


+  $(-A  sin  3)$(A  cos  &-Jt) 


+  $(-A  sin  (tt/3  -  3))$(A  cos  (tt/3  -  3)  -</t). 


Bartholomew  (1961)  took  the  partial  derivative  with  respect  to  3  and 
noted  that  for  a  fixed  value  of  A,  the  power  function,  which  is  periodic 
with  period  2tt,  is  increasing  for  3  €  [-5tt/6,  tt/6]  and  is  decreasing 
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for  3  €  [tt/6,  7tt/ 5] .  Thus,  it  has  a  maximum  at  3  =  tt/6  (the  middle  of 
Hj),  a  minimum  at  3  =  7n/6  (the  middle  of  H*)  and  since  it  is  sym¬ 
metric  about  3  =  tt/6,  it  has  two  equal  minima  within  H-j  at  3  =  0  and 
3  =  tt/3. 

The  partial  derivative  of  (3.1)  with  respect  to  A,  evaluated  at 
A  =  0,  is  00/t) (v/3/2  +  yt/ (2tt) )  cos(3-tt/6),  which  is  positive  for 
3  €  (-tt/3,  2rr/3)  and  negative  for  3  6  (2tt/3,  5tt/3).  This  might  lead  one 
to  conjecture  that  the  power  of  T01  15  increasing  in  A  s  0  for  fixed 

3  €  (-tt/3,  2r r/3).  We  have  not  been  able  to  establish  that  using  (3.1). 
However,  the  results  given  in  this  section  (Corollary  3.2)  imply  that  it 
is  increasing  in  A  i  0  for  fixed  g  6  (-tt/6,  tt/2).  We  have  also  shown 
that  this  is  the  case  for  k  =  3  and  3  €  (-tt/3,  2tt/3)  using  techniques 
similar  to  those  employed  in  Section  2,  but  because  of  their  special 
nature  these  arguments  are  not  given  here.  Applying  the  results  of  Theorem 
3.3  and  those  concerning  the  sign  of  the  derivative,  at  A  =  0,  of  the 
power,  we  know  that  for  3  6  (-tt/2,  -tt/6)  U  (2rr/3,  5tt/6)  the  power  func¬ 
tion  first  decreases  in  A  and  then  approaches  1  as  A  ->  00 . 

Bartholomew  (1961)  also  considered  the  case  k  =  4,  but  the  expres¬ 
sion  for  the  power  function  is  quite  complicated.  Several  values  of  the 
power  function  for  k  =  4  are  computed  there  and  results  like  those 
obtained  for  k  =  3  are  conjectured  to  hold  for  k  =  4,  but  no  further 
analysis  of  the  power  function  is  given  for  k  =  4. 

Remark.  If  a  ^  -1  and  v  f  Rk  then  the  distance  from  v+aEw(v[  ) 
to  is  the  same  as  the  distance  from  v  to  H-j . 


Proof.  Using  B.l  and  B.2  on  page  131  of  Barlow  et  al .  (1972),  it  can 


Theorem  3.1.  Let  v  €  R  .  As  a  function  of  a  €  (-«,*), 

(v+  a£^(v  (  H^))  is  nondecreasing  and  ttq-|  (v+ aE^(v  j  H*^))  is  non¬ 
increasing. 

Proof.  Write  v  + aE^(v  J  H.j )  =  v  -  E^(v  |  H-| )  +  (a+1  )E^(v  |  H  ) ,  set 
=  v-  Ew(v  |  )  =  Ew(v  |  H*W)  and  ( a  =  Ew(v  |  H] )  and  note  that  by  (2.1), 

(nQ,la)w  =  0.  Set  S  =  HjW  =  (x  6  Rk:  rf=1  w.x.  =  0}  and  C  =  CQ1  =  H^S. 

i  W  i  u  .  C  v 

Since  S  =  Hg,  Cg-j  ©  S  =  .  Applying  Theorem  2.10,  W(AW) 

-  Ppi  +6|j  w^oi  - is  nondecreasing  in  6.  Thus,  (n0+6^i)  is  non¬ 
decreasing  in  5^0.  For  6  s  0,  consider  (ng+(-6)  (-u) ) ,  which  is 
nonincreasing  in  -6  if  -pi  €  C*^  =  (H-j  H  HgW)*W  =  H*W  ©  Hg.  Barlow  et  al . 
(1972,  p.  49)  show  that 

(3.2)  H*W  =  (x  €  Rk  :  I1.  ,  w.x.  i  0  for  i  »  1,2,  —  ,k-l 


Jj=l  J  J 


and  Z.a]  w.x.  =  03. 


Since  (-H] )  n  HjW  c  H*W,  -pi  +  (£k=1  wjUj)ek  €  H*W  and  -u  €  H*W  ©  HQ, 
the  first  claim  is  established. 

For  the  second  conclusion,  write  v+aE^(v|  H*W)  as  pig  +  (a+1  )pi  with 

Ug  =  Ew(v  j  )  and  pi  =  Ew(v  i  H*W)  and  note  that  (nQ,pi)w  =  0.  Let  S 

and  C  be  as  in  the  first  part  of  the  proof.  Now  n  6  c  Cg^  =  (H^flS) 

*W 

=  H.j  ©  Hg,  and  applying  Theorem  2.10,  ttq-|  (pig+6n)  is  nonincreasing  in 

*W 

6  i  0.  For  6^0,  consider  nQ]  (pig+  ( -6 )  (-pi) )  and  note  that  -pi  €  -H-j 

.  u  ..... 

c -C  .  Thus,  applying  Theorem  2.10  again,  rrg^(pig+6u)  is  nomncreasing 

for  ft  s:  0.  The  proof  is  completed. 


Theorem  3.1  can  also  be  established  using  the  results  in  Robertson  and 
Wright  (1982).  They  considered  two  relations  on  R  defined  by  x<y 


*w 

provided  y-x  €  and  x  «  y  provided  y-x  €  -H^  and  proved  that  if 
either  pi  <  | pi'  or  pi  «  pi'  then  tt^  (pi)  s:  tt01  (pi' )  (i.e.,  nQ1  ( • )  is 
isotone  with  respect  to  both  £  and  «).  For  a-|  ^  c^.  v+a-|E^(v|  ) 

<v+a2E^(v|  H.| )  and  v+o^E^vj  H*^)  «  v+a-|E^(v(  H*^)  and  the  results 
of  Theorem  3.2  follow  immediately.  However,  this  approach  does  not  yield 
the  analogous  results  for  T^- 

Corollary  3.2.  If  pi  6  -H*W  ©  Hq,  then  ttq^(6u)  is  nondecreasing  in 
6t(-°°,®).  Furthermore,  is  unbiased. 

Proof.  By  hypothesis,  pi  =  v+v'  with  v  €  -H*W  and  v'  €  Hq. 
Examining  the  definition  of  ,  it  is  clear  that  ttq-|(6u)  =  ttq.j(6v). 
Applying  the  second  conclusion  of  Theorem  3.1  with  v  replaced  by  -v, 
we  see  that  tTq^C&v)  is  nondecreasing  in  6  € 

.  u 

In  the  proof  of  Theorem  3.1,  we  saw  that  c:  -H^  ©  Hq.  So  for 
pi  €  H-j-Hq,  ttq1(pi)  2:  nQ1(0-u).  The  proof  is  completed. 

*W 

Remark.  Since  H-|  ©  Hg,  Tiffin)  is  nondecreasing  in  6  6  (-®,® 

if  pi  €  H, . 

The  limiting  behavior  of  ,  in  the  directions  considered  in  Theorem 
3.1  and  Corollary  3.2,  will  be  discussed  next.  A  slight  generalization  of 

Theorem  3.4  of  Barlow  et  al.  (1972)  will  be  used  to  obtain  the  limits  of 

.  •  .  k 

interest.  Following  their  notation,  we  define  for  pi  6  R  , 

A2  =  a2(u)  =  w.^-pT)2  with  u  =  wini/L^=1  w.. 

We  first  note  that  the  statement  of  the  theorem  given  there  needs  some  clar¬ 
ification.  Consider  the  following  example:  Let  w  =  e^,  let  v  6 


with  Uv-vii  >  0  and  let  nn  =  (-1)  nv.  Then  A  (un)  =  n  A^(v)  ->  «>  and 
(u  .-u  )/ (u  .-u  )  =  (v--v)/(v.-v)  for  1  <:  i,  j  «:  k.  However,  for  odd 
n  the  number  of  distinct  values  in  E(nn|H^)  is  1,  i.e.,  E(nn|H^) 
is  constant,  but  for  n  even,  E(un(  )  =  n-v  which  is  not  constant 

since  Hv-v||  >  0.  Hence,  the  hypotheses  of  their  Theorem  3.4  may  hold,  but 

may  not  be  constant.  By  making  only  slight  modifications  of  their 
proof,  one  can  prove  the  following  generalization. 

i/ 

Theorem  3.3.  Let  v,0,rin  €  R  with  rin  0  as  n  ®.  If  wn  =  apv 
and  un-iln  =  by2(e-0+rin)  with  anbn  ->  ®  as  n  ->  »,  then  ^(i^)  1 

provided  0  /£  H*W  ©  Hg,  and  ^q]  (l^n )  0  if  8  €  (H*W©Hg)°,  where  A0 

denotes  the  interior  of  A. 

Barlow  et  al .  (1972)  applied  this  result  to  show  that  T^  is  consis¬ 
tent  for  n  £  H*W  ©  Hg .  It  also  gives  the  radial  limits  of  the  power  func¬ 
tion  in  certain  directions.  These  limits  are  obtained  by  setting  v  =  0 
in  the  following: 

k  *u 

Corollary  3.4.  Let  ta,v  €  R*.  If  u  £  H^  ©  Hg,  then 
11m  TToi  (v+6ja)  =  1-  If  n  €  (H*W©Hg)°,  then  1  im^  nQ1  (v+6u)  =  0. 

Proof.  The  result  follows  from  Theorem  3.3  by  setting  v  =  w,  0  =  u, 

"n  *  (v-v)/6n,  an«l,  bn  -  6*  with  6„ 

*w 

One  can  obtain  the  form  of  1  im^  TTg^  (v+6ta)  with  n  €  9(H^  ©Hg), 
but  such  limits  play  no  central  role  in  this  work  and  are  tedious  to  develop, 
so  they  are  not  included.  Next,  we  study  the  limits  of  the  power  function 
along  lines  parallel  to  H^  and  H*^,  that  is,  we  consider 

11ma>±-TT01(v4aEW(vlHl,)  and  limc^-TT01(v+aEW(v*HlW))-  Directions 


V  V 


V  ’ 


- .  - .  . 
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parallel  to  H-j  will  be  discussed  first. 


Corollary  3.5.  Let  v  €  R*.  If  v  €  H*  ©  Hq,  then 


ttqi  (v+o£w(v  i  H.]))  =  nQ1(v).  If  v  £  H*  ©  HQ,  then 

lima>~  ^oi(v+aEW(v  !  H))  =  1  and  11 V—  n01(v+aEW(vl  Hl})  =  °* 


Proof.  Using  (3.2),  we  characterize  H*^  ©  Hq  as  follows: 


(3.3)  H*W  SHQ  =  {x  €  Rk  :  ll=1  w.x.  *  (E’=1  w^x  for  i  =  1 ,2,  • • • ,k-l ] . 


Using  the  minimum  lower  sets  algorithm  for  computing  E(x  [  )  (cf.  Barlow 


et  al.  (1972,  p.  76)),  we  see  that  Eu(x  (  H, )  €  H  if  and  only  if 

W  1  o 


J*W  /r\  u  Ti  .  .  r  LI*W 


x  €  H*  ©  Hg.  If  V  €  H*  ©  Hq,  then  the  first  conclusion  follows  from 


★w 

the  fact  that  tt^  is  invariant  under  constant  shifts.  If  v  £  ©  Hq, 

then  Ew(v(H1))6H0  and  since  H,  n  (H*W  ©  HQ)  =  HQ,  Ew(v  |  H] )  £  H*W  ©  HQ. 

Applying  Corollary  3.4  gives  the  second  conclusion. 

For  the  last  conclusion,  we  consider  tt01 (v+(-a) (-Ew(v  |  H-j))),  The 

desired  result  follows  from  Corollary  3.4  by  showing  that  -Eyj(v  |  H^)  is 

u 

in  the  interior  of  H*  ©  HQ,  which  is  characterized  by  making  the  ine¬ 
qualities  in  (3.3)  strict.  Now  -Ew(v  |  H^ )  £  HQ  and  has  nonincreasing 


coordinates,  and  for  any  x  f  T  with  these  two  properties. 


Z1.  ,  w.x./Z1.  ,  w.  is  nonincreasing  in  i  and  equals  x  for  i  =  k 
J"  '  J  J  J  ^  J 

vk"l  ..  „  ^  ...  ^  *7?  w  T 


Furthermore,  w.x./E."  w.  >  x,  for  if  not  x  is  constant.  The 

J  *  3  0  J  I  J 

proof  is  completed. 


*W  *W 

We  now  consider  limits  along  lines  parallel  to  H^  .  If  v  €  H-j  , 


v  +  aEw(v[  H*  j  =  (a+l)v  and  yields  radial  limits  as  a 


v  If  H|W,  then  Ew(v  |  H*W)  €  3(H*W  ©  HQ)  and  1  im^  ttq1  (v+aEw(v  i  H*W) ) 


is  the  type  of  limit  discussed  after  the  proof  of  Corollary  3.4. 


Corollary  3.6.  Let  v  €  R^.  If  v  €  ,  then  tt^  (v+aE^(v  |  H*^)) 

=  nQ1(v).  If  v  it  then  lim^^  ttq1  (v+aEu(v  |  H*W))  =  1. 

Proof.  Since  E^(v  {  H*^)  =  v-E^(v  |  H^),  we  see  that  E^(v  |  H*W)  =  0 
if  and  only  if  v  €  .  The  first  conclusion  is  now  clear. 

For  the  second  part,  assume  v  It  H^.  Using  (3.2)  and  (3.3),  we  see 
that  x  €  H*W  ©  Hg  <=>  x-x  €  H*W.  The  proof  is  completed  by  applying 
Corollary  3.4  provided  -Ew(v  {  H*W)  £  H*W  ©  HQ.  If  -Ew(v  |  H*W)  6  H*W  ©  HQ, 

then  -Ew(v|  H*W)  €  H*W.  But  (-H*W)  n  H*W  =  {0},  and  if  Ew(v  !  H*W)  =  0, 

then  v  €  H-j . 

We  now  turn  our  attention  to  the  study  of  tt^„»  The  power  function  of 

T-|2*  For  k  =  2  rejecting  in  favor  of  for  large  values 

of  T^  is  equivalent  to  rejecting  if  X^-Xg  is  large.  This  test  is 
known  to  be  unbiased,  UMP  and  to  have  a  power  function  which  is  nondecreas¬ 
ing  in  For  k  =  3  and  w^  =  w2  =  w^,  one  can  employ  the  same 

techniques  used  by  Bartholomew  (1961)  to  show  that 

exp(-^-A2)  b-tt/3 

(3.4)  TT,p(n)  =  P  [T ■» o >  t]  = - - - f  *(A  sin  0,t)d0 

u  ^  u  2tt  0-tt 

+  $(-A  sin  0 cos  0) 

+  §(-A  cos(B-hr/6)-y/t)5(A  sin(0+n/6)), 

with  A  and  0  defined  as  before.  It  is  not  difficult  to  show  that,  for 

fixed  A,  tti 2 ( * )  1S  antisymmetric  about  0  =  rr/6  and  0  =  7n/6.  Based 

on  Bartholomew's  work  on  with  k  =  3  and  w-j  =  w?  =  W3»  one  would 

conjecture  that  is,  for  fixed  A,  decreasing  for  0  €  [-5rr/6,  tt/6] 

and  increasing  for  0  €  [tt/6,  7tt/6].  We  have  not  established  this 


?  '  B  »  t »  l1."  IT’  r-v*  V  ■.  ■  ■  ■  «j  mi  y '  -  I  ^  n;viV.  ■,-»  v 

•  f  »  *  .«  _•  _ *  „►  “  .*  «  ^  .  •  _  *  .  "  «  “  » 
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analytically,  but  have  numerically  obtained  the  value  of  for 

several  values  of  0,  A  and  t.  This  partial  derivative  appears  to  be 
negative  on  (-5tt/6,  tt/6)  and  positive  on  (tt/6,  7tt/6). 

The  partial  derivative  of  n^*  with  respect  to  A,  evaluated  at 
A  =  0  is  -0(v/t) (3t/(2n) ) ^  +  1/2)cos(0-tt/6) ,  which  is  negative  for  0€ 
(-tt/3,  2tt/3)  and  positive  for  0  €  (2n/3,  5tt/3).  As  one  might  expect, 
this  behavior  is  opposite  to  that  of  . 

For  arbitrary  k,  we  apply  the  results  in  Section  2  in  our  study 
of  tt^. 

Theorem  3.7.  Let  v  6  Rk.  As  a  function  of  a,  tt^2(v+ a£w(v  {  H^)) 
is  nonincreasing  for  -«  <  a  <  ®  and  tt-|2(v+ aE^(v  [  H*W))  is  nondecreas¬ 
ing  for  a  £  -1 . 

Proof.  Write  v+aEw(v|  H^)  =  Ew(v  [  H*W)  +  (a+l)Ew(v|  H-j )  =  uQ  +  (a+1  )u 

with  (u0,u)w  =  0.  Set  S  =  HjW  =  [x  €  Rk :  Ek=1  w.x.  =  0)  and 

*W  +W 

C  =  C-j 2  =  H-j  CS.  Applying  Theorem  2.10,  we  see  that  for  n  e  =  C-^* 

r 

P6u  W<V  =  n12^0+^  1S  nonincreasing  in  6  for  6^0.  For  6^0, 
consider  which  is  nondecreasing  in  -6  since 

£■  u  -  r*W 
-u  €  -H-j  -  -C-|2‘ 

For  the  second  conclusion,  v+aEw(v[  H*W)  =  Ew(v  {  )  +  (a+l)Ew(v[  H*W) 

★w 

=  Uq  +  (a+l)w  with  (n0»u)w  -  0.  Since  u  €  H-j  =  we  apply  Theorem 

2.10  to  show  that  n^luQ+Su)  is  nondecreasing  in  6  s  0.  The  proof  of 
the  Theorem  is  completed. 

Comparing  Theorems  3.1  and  3.7,  we  see  that  ttq.| (v+ aE^(v  f  H*W))  is 
monotone  for  <  a  <  00 ,  but  tt^ 2(v+ aE^(v  (  H*^))  is  only  claimed  to  be 


monotone  for  a  ^  -1.  We  consider  an  example  to  show  that  the  second  con 
elusion  of  Theorem  3.7  is  not  valid  for  -®  <  a  <  ®. 

Example.  Let  k.  =  3,  w  =  e^,  v  =  C/2/ 2, --72/2,0)  €  H*.  (Recall,  H 
is  characterized  in  (3.2).)  Now  v  +  aE(vi  H*)  =  (a+l)v,  and  we  will  show 
that  n12((a+l)v)  is  not  monotone  in  (-®,-l).  If  it  were,  then 
tt12(6(-v))  would  be  monotone  for  6  >  0.  However,  the  g  corresponding 
to  -v  is  tt/2  and  for  such  g,  <  0-  Hence,  the  power 

decreases  for  6  small  and  positive,  but  applying  Corollary  3.11,  we  see 
that  !  im^  rr12(&(-v))  =  1 . 

The  above  example  is  interesting  for  several  other  reasons,  also. 

(1)  For  the  v  chosen,  v+aEw(vf  H*)  =  (a+1  )v  and  so  we  see  that 
tt.j2(6v)  is  not  monotone  for  6  <  0  (see  the  next  corollary). 

(2)  It  shows  that  is  not  antitone  with  respect  to  the  partial  order 

«,  discussed  in  Robertson  and  Wright  (1982),  that  is,  if  0  s  6<  6 
then  6(-v)  «  6' (-v)  but  n^(-6v)  may  be  less  than  n^2(-fiyv). 

(3)  It  shows  that  T^  is  biased.  Along  the  ray  [6(-v)  :  6  £  0],  the 
power  decreases  for  small,  positive  6  and  so  the  level  of  signifi¬ 
cance  is  at  least  (and,  in  fact,  is  equal  to)  rr-j2(0)  =  tt^2(0(-v)) 
>tt-|2(6(-v))  for  some  6  >  0.  We  will  consider  the  question  of  the 
unbiasedness  of  T^2  in  more  detail  later  in  this  section. 

Corollary  3.8.  If  n  €  H1 ,  then  n12(6fa)  is  nonincreasing  for 
ft  €  (-®,®).  If  n  6  H*W  ©  Hq,  then  tt12(6m)  is  nondecreasing  for 
6  €  [0,®). 

Proof.  Corollary  3.8  follows  from  Theorem  3.7  just  as  Corollary  3.2 
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follows  from  Theorem  3.1. 

To  obtain  the  limits  of  in  the  directions  considered  in  Theorem 

3.7  and  Corollary  3.8,  we  establish  the  analogue  of  Theorem  3.3  for  tt 

Theorem  3.9.  Let  v.e.r^  €  Rk  with  nn  ->  0  as  n  -»  «.  If 

Wn  =  anv  and  un'wn  =  b,!/2(9"^+,rln)  w'th  anbn  ®  as  n  00 .  then 

n12^n^  1  Provided  6  jf  H]  and  Tr12(nn)  ->  0  if  6  6  hJ. 

Proof.  The  proof  of  the  first  conclusion  is  very  similar  to  the  proof 
given  in  Barlow  et  al.  (1972)  for  their  Theorem  3.4.  The  LRT  of  versus 
H2  rejects  if  T^2  >  t.  In  this  case,  T ^ 2/ ( bn )  converges  in  prob¬ 
ability  to  Zk=1  vi  (E(e  1  H^-9  )2,  which  is  positive  if  0  £  H  .  Of  course, 
t/(anbn)  ->  0  and  the  first  conclusion  is  established. 

For  the  second  conclusion,  tt^)  =  Pp  ^{x  :  JjEw(x  |  >  t). 

But  E^(*  |  H-j)  =  Ey('  J  H-|)  where  V  is  a  kxk  diagonal  matrix  with 
Vii  =  vi  ^  i  =  1  *2, •  •  •  ,k,  |lEv(x!  -  I|Ev(a’/2x  (  H,  )-a’/2x|lg. 

and  if  x -t^.W"1)  then  aj/2x  ~  nva^^.V"1 ) .  So,  n12(Up) 

P0,V*-X  :  Hi  )“x-ap/,2Wn(lv  >  t}.  By  the  hypotheses  of  the 

theorem,  aj/2^  =  (anbn)1/2(e-3+Tin),  aRbn  ®,  ->  0  and  e-ff  has 

strictly  increasing  coordinates.  Thus,  for  each  x  €  Rk,  there  exists  an 
n(x),  with  x+ay2un€H1  for  all  n  st  n(x).  Hence, 

HEv(x+ay2un  i  Hi }  -  x  -  ap/2unliv  =  0  forall  n  s  n(x) , 
and  because  t  >  0,  the  desired  result  is  established. 

Corollary  3.10.  If  n..  =  ny..  with  y^  >  0  for  i  =  l,2,***,k,  then 
T-j2  is  consistent  for  all  y  |f 


Proof.  This  result  follows  from  Theorem  3.9  by  setting 
2  2 

V  =  (Y^.  "»Yk/<*k)»  e  *  y,  r\n  E  0,  bn  e  1,  and  an  =  n. 

Corollary  3.11 .  Let  y,v  €  Rk.  If  y  £  ,  then  lim6_^  n^v+fiy) 

=  1,  and  if  y  €  H^,  then  1  imfi_^  ^(v+bn)  =  0. 

Corollary  3.11  follows  immediately  from  Theorem  3.9  and  with  v  =  0, 
gives  the  values  of  certain  radial  limits.  Because  the  radial  limits  of 
for  y  €  dH.j  are  of  interest  in  our  study  of  the  bias  of  T12’  we 
need  to  obtain  the  value  of  these  limits.  For  y  €  H.j ,  let 

1  ^  <  d2  <  ' '  ’  <  Jh  =  k  be  definec!  n-j  =  *  *  •  =  Uj  <  =  ’  ’  ’ 


,+l  . V  ^ 


(3.5)  C/  (y)  =  {x  €  R  :  x.  £  x„s---  *:  x.  ,x.  +1^  •••  £  x 

•  a  J  i 


* • *  »X^  .  S  •  • .  £  X  .  ) 

-h-1  1  Jh 

and  set  G-j  -  £l  »2»*  *  *  ,j-|3»  G2  ~  { d-j^l » ■  ’  ’  >  ’  ’  *  >G^  -  »j^3» 

The  G  are  the  level  sets  of  y, 

JL 


Theorem  3.12.  Let  v  f  R  and  u  6  H^.  Then, 


(3.6)  lim6^  tt12(v+6u)  =  Pv>wC|iEw(x|  C'(y))-xjjJ  >  t]. 

Proof.  If  y  €  H^,  then  C '  (y)  =  Rk,  the  r.h.s.  of  (3.6)  is  zero 
and  (3.6)  follows  from  Corollary  3.11.  Suppose  y  6  and  consider 
Eu(x+v+6y!H1).  For  a  fixed  x,  x.  +  v.  +  6y .  -  x .  -  v .  -  6y .  ->  ®  as  6  ->  ® 
for  i  €  G^,  j  6  G^  with  l*  <  A.  So  for  each  fixed  x,  there  exists  a 
6(x)  with 


V  V  V  V  / 


.  •  •  ' .-w ,  .  • .  •  .v.s  ;. 
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maxi-€G1(xi+Veui)  <  mini€G2(xf+vi+6w1)  S  maxi€G2(x1+vi+6^i) 


<  •••  <  min.  G  (x.+v^^.) 
h 


for  6  £  6(x).  It  follows  from  the  minimum  lower  sets  algorithm  (cf.  Barlow 

et  al.  (1972,  p.  76))  that  for  in  €  G.,  min.,r  y.  s;  y .  £  max.^  y.  and 

0  1  l€GX  1  70  l€Gx  1 

that  since  6u  is  constant  on  the  level  sets  of  u,  Ew(y  +  6u)C' (pi)) 

=  Ew(y)C'  (u))  +  6u.  Since  C'(u)^H^,  Ew(x  +  v+  6pt  |  C'  (u) )  =  Ew(x  +  v+  6pi  |  ) 

for  6  £  6(x).  Hence,  for  each  x,  jjE^(x  +  v+  i  H-j )  -  x  -  v  -  6u|Jw 
->  !IEw(x|  C#(u))-x|!w.  Thus, 


tt12(v  +  6u)  =  P0jW[x  :  |]Ew(x  +  v  +  6h  1  H-j)  -  x-  v-6ujj^  >  t]  -> 
P0,^x:  !IE«(X+VI  C'  (ja))  -  x  -  v||y  >  tj 
-  Pv>„[!|E„(x|c'W)-x^>t]. 


The  proof  is  completed. 


By  taking  v  =  0  in  Theorem  3.12,  we  obtain  the  radial  limits  for 

pi  €  SH^.  Corollary  2.6  of  Robertson  and  Wegman  shows  that  the  r.h.s.  of 

2 

(3.6)  with  v  =  0  is  a  weighted  sum  of  \  tail  probabilities,  and  the 
remark  on  p.  148  of  Barlow  et  al.  (1972)  shows  that  the  weighting  constants, 
i.e.,  level  probabilities,  in  this  case,  are  convolutions  of  those  for  a 
total  order.  We  will  compute  some  values  for  this  limit  when  we  study 
bias.  However,  since  [!Ew(x  1  Hj )  -  x||w  is  [|Ew(x  |  C'  (u))  -  x[J.w 

and  so,  Pv>wCi!Ew(x  |  C'(u))  -  x||w  >  t]  *  n]2(v). 


Next,  we  study  the  limits  of  the  power  function  along  lines  parallel 


to  H1  and  H*W.  Directions  parallel  to  H1  will  be  discussed  first. 

If  v  €  H] ,  then  v+aEw(v(  H] )  =  (a+l)v,  which  yields  a  radial  limit  as 
a  ->  ®.  So  we  may  suppose  v  if  H^. 

Corollary  3.13.  If  v  £  ,  then 

1  ima-*=°  tt-|2^v+  aMv  !  Hl^  =  Pv,W^MX  l  (EW(v  ^  Hi  ^ )  )“xliw  > 

If  v  €  H*W  ©  HQ,  then  n] 2(v  +  aEw(v  |  H] ) )  =  n12(v) ,  and  if  v  £  H*W  ©  HQ, 
then 

^V-co  n12(v+aEw(v|  H] ) )  =  1. 

Proof.  The  first  part  of  the  result  follows  from  Theorem  3.12.  For 

*W 

the  second  part,  we  recall  that  Ey(v  (  H-j)  €  Hq  if  and  only  if  v  €  ©Hq. 

u  *W 

So,  if  v  €  H*w  e  H0,  then  rr12(v+  aEw(v  {  H] ))  =  rr12(v).  If  v  £  H]  e  HQ, 

then  -Ew(vi  H-j )  £  H]  because  (-H^  D  H]  =  HQ.  Applying  Corollary  3.11 

gives  the  desired  conclusion,  and  the  proof  is  completed. 

Directions  parallel  to  H*W  are  considered  next.  We  may  assume 
v  ?  H*W,  and  in  fact,  since  tt-j 2  is  invariant  under  constant  shifts,  we 
may  assume  v  £  H*  ©  Hq. 

Corollary  3.14.  If  v  €  H] ,  then  n12(v  +  aEw(v  [  H*W) )  =  tt12(v) .  If 
v  £  H^ ,  then  lirn^  Tr12(v  +  aEw(v  i  H*W))  =  1.  If  v  £  H-j  U  (H*W©HQ), 
then  1  im^.^  TT12(v  +  aEw(v  |  H*W) )  =  1 

Proof.  If  v  6  H-j ,  then  Ew(v  (  H*W)  =  0  and  so  the  first  conclu¬ 
sion  is  immediate.  If  v  £  H-|,  or  equivalently  E^(v  i  H^  )  £  H-j,  then 

*W\ 

appealing  to  Corollary  3.11,  1  ttt2^v+  aEW^v  ‘  ^1  '  =  "^e  ^ast 
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conclusion  also  follows  from  Corollary  3.11,  if  we  can  show  that 


-Ew(v|  H*W)  £  H1 .  Suppose  n  =  -Ew(v  [  H*W)  6  H] .  Then  L^=1  w^  =  0,  t}. 
is  nondecreasing  and  ri  £  HQ  (Ew(v  (  H*W)  6  HQ  <=>  Ew(v  1  H*W)  =  0  <=*>  v  6  H-j), 
So,  WjT]j  is  nondecreasing  in  i  and  equal  to  zero  for  i  =  k.  Since 

\  t  o,  e!=1  Wjtij  <0  for  i  =  1 ,2,* •  •  ,k-l .  If  G1  =  £l ,2, *--,0^3  is 
the  first  level  set  for  Ew(v  i  H-j ) ,  then  31  <  k  and  w.t^ 

■  Ej'l  Wj(EH(vl  HlW  ■  Zj=1  VA'l  W0V0  =  °-  ™S 


contradiction  completes  the  proof. 


We  have  already  noted  that  is  biased  and  we  now  wish  to  examine 


the  amount  of  bias.  In  the  case  k  =  3  with  w  =  e^,  the  level  of  sig¬ 


nificance  is  tt^(0)  =  t]/3+  P[x^ >  t]/2.  Partition  R3  into  four 


sets  depending  on  the  number  of  level  sets  in  the  projection  onto  H. 
Specifically,  with  x*  =  Ew(x  (  ) ,  let 


C,  =  £x:  x*<x^<x*3  (=  Hp ,  C„  =  {x  :  x*=x*<x*} 


12  3 


4>  4>  ■ 

•  —  v  ✓  V  A  I 


1  2  3- 


C3  =  {x  :  x*<  x*  =  x*J  and  =  {x  :  x*  =  x£  =  xp  (=  H*  ©  HQ) 


We  have  seen  that  inf^  tr-^n)  =  0  (cf.  Corollary  3.11)  and  that 


inf|a€C  =  (cf.  Corollary  3.8).  It  will  be  shown  (cf.  Theorem 


3.15)  that  inf  q.  tt12(u)  =  inf  ^  tt-^Gj)  =  P[x^>t]/2,  and  so  by  the 
2  3 


rlnr  2 . 


continuity  of  n^,  inf^^  tt^2(u)  =  2  PCx^t].  In  the  case  being 


sidered,  the  5%  critical  value  for  T1?  is  4.578  and  P[xf^  4.578]/2 


=  .0162,  which  gives  some  idea  of  the  amount  of  bias.  (Larger  k  will  be 


discuss'd  later.) 


N»',v 


Returning  to  the  case  of  arbitrary  k,  we  partition  R  into  M  =  2 
subsets  depending  on  the  level  sets  of  x*.  Let  =  {x  :  x*< x£ <  • • • < x£] 
(=  H°),  C2  =  £x:  x*  =  x*  <  x*<  --^x*},  C3  =  (x:  x*<x*=  x*<xj<  *..<xj 
•••»  CM  =  [x  :  x*  =  •••  =  xp  (=  H*W  ©  Hq).  For  x  €  C^,  let  C'(x*)  be 
defined  as  in  (3.5)  and  note  that  C'(x*)  is  the  same  cone  for  each  x  €  C 
Set  =  C' (x*)  for  x  €  C^. 

Theorem  3.15.  With  C.  and  C'.  defined  as  above, 

(3.7)  ^12^)  “  ^  ^ 

Proof.  We  consider  =  H*W  ©  HQ  first  and  note  that  C'M  =  H-| . 
Equation  (3.7)  follows  from  Corollary  3.8.  Fix  i  <  M  and  let  n  €  . 

Now,  as  in  the  proof  of  the  Remark  preceding  Theorem  3.1,  |j+aEw(u|  H1 ) 
has  the  same  level  sets  as  pi  for  a*  -1.  Applying  Theorem  3.7, 
n12(u  +  aEw(n  (  H1 ))  ^  rr12(n)  for  a  0  so  that 


i'Vq  TT12(|j)  =  1nfu€Cj11n,a-I,]2(M+aEW(wI  V1 
=  inVi  pn,w[l|Ew(xlci)'*l:w>t] 

by  Corollary  3.13. 

Let  =  (^  +1 , • • • )  and  W^  be  the  <3x"3x_i ) x  Ux-Jx  ^ 

X™  1  X  "1 

diagonal  matrix  with  diagonal  elements  w.  ,,•••, w.  for  X=l,2,***,h. 

JX-1  °X 

Now,  jjEw(x  [  Cp-xJly  =  ZJJ=1!|EW  (x£|  H]  jX)-xJJ^  ,  which  is  a  sum  if  inde- 

X  5  X 

k  th 

pendent  random  variables  on  R  .  The  distribution  of  the  X  summand 


could  be  thought  of  as  indexed  by  (u,W)  or  (ux>Wx).  In  the  latter  case, 
we  can  apply  (3.7)  with  i  =  M,  to  see  that 


P^,w/llEWx(xX  5  H!  u)-XJ>Wx  >  t]  *  P0,w/llEw/XX  l  H1  ,*}-XAx  >  «  Since 

Ew(n^  [  H-j  j)  is  constant.  Under  both  probabilities  P  w  and  Pq  w, 

'  '  2  th 

j|Ew(x  j  C^)-x||^  is  a  sum  of  h  independent  random  variables  with  the  JL — 

summand  stochastically  larger  under  P  w  then  under  PQ  l  -  l,2,-*-,h 
2  ^  * 

So  I|Ew(x  |  C'- )  -x[jw  is  stochastically  larger  under  P^  w  (cf.  Proposition 
C.l,  p.  485,  Marshall  and  01  kin  (1979)).  Hence, 

infu€Ci  Pn,W^‘EW^X  ^  ci)_x!‘w  >  ^  *  P0,W^‘Mx  f  ci}-xliw  > 

and  the  reverse  inequality  follows  from  the  fact  that  C.  is  a  cone  and 
P  u  is  continuous  in  |i. 

Corollary  3.16.  inf^H  tt12(h)  = 

Proof.  By  the  continuity  of  inf^  tt12(u)  =  inf^  tt12(h). 

Fix  i  >  1,  then  there  is  some  j  with  xt  =  x^  for  all  x  f  C.  and 

Cj+1  =  [x:  x*<  <  x*  =  x*+]  <  ...  <  x*j.  Hence,  C'  c  C'+1 

=  [x  €  Rk:  x.  z  xj+1},  ||Ew(x|  C'^-x^  *  ||Ew(x  |  Cj+1)-x[|w  for  all  x  6  Rk 

anC^  Po,WtllEW(x  i  ^ _XH w  >  ^  ^  p0,WniEW(x  I  ^j+l^-xl*w  >  ^°> 

1pWhi  n12^  =  infl^jsk  P0,W^‘EW^x^  Cj+l^"xIIw  >  But’  EW(xE  Cj+l^i 
=  xi  for  i  f  j,j+l  and  (Ew(x  !  C'.^)j,  Ew(x  j  C'-+1)j+1)  is  the  projec- 

2 

tion  of  (xj»xj+])  onto  (y  €  R  :yi  s  y23  with  norm  defined  by 

||y||2  =  vry^  + Wj+.jy|.  Using  Corollary  4.2  of  Robertson  and  Wegman  (1978) 

and  the  fact  that  for  a  total  order  and  any  weights  P ( 1 , 2 )  =  P ( 2 , 2 )  =  1/2, 
we  see  that  p0>wCJ{Ew(x  I  C'.+])-x|jw  >  t]  =  P[x^  >  t]/2.  The  proof  is  com¬ 
pleted. 
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It  is  of  interest  to  examine  the  "amount"  of  bias  in  as  k 

increases.  So  with  t  the  5%  critical  value  for  with  w  =  e^  and 

k  =  3, 4, 5, 6,  we  computed  the  infimum  in  Corollary  3.16.  As  was  noted 
before,  for  k  =  3  the  infimum  is  .01620,  for  k  =  4  it  is  .00648, 
for  k  =  5  it  is  .00281  and  for  k  =  6  it  is  .00128.  (The  critical 
values  are  taken  from  Robertson  and  Wegman  (1978)).  The  infimum  is  approx¬ 
imated  by  n^(n)  with  n  at  a  large  distance  from  Hq,  but  close  to  . 
For  practical  purposes  it  is  also  of  interest  to  compute  tt^  at  u  near 
H^,  but  at  a  "reasonable"  distance  from  Hq.  We  first  consider  k  =  3 

and  w  =  e3,  for  in  this  case  the  powers  can  be  obtained  numerically. 

Because  of  the  proof  of  Corollary  3.16,  we  will  compute  +  |  H^)) 

for  various  a  and  v  chosen  so  that  E(v|  H-j)  has  one  level  set  with 
two  elements  and  the  other  has  one  element.  Table  1  gives  the  values  of 
tt12(v.  +aE(vi  \  H.j))  for  vi  =  p^/AU-i-j)  with  =  (2,1,2)  and 
H2  =  (1.5, 1,2)  and  t  =  4.578,  the  5 %  critical  value. 

We  observe  that  the  bias  of  T^»  even  f°r  k  =  3,  is  large  enough 
to  be  of  practical  significance.  For  k  =  5,  w  =  e^,  t  =  7.665  (the  5% 

critical  value  of  T^)  and  v  =  piM(u)  with  |j  =  (3, 1,3, 4, 5), 

Tr12(v  +  aE(v  I  H^))  are  estimated  by  Monte  Carlo  techniques  with  10,000 
repl ications.  These  values  are  given  in  Table  2.  We  notice  that  the  bias 
is  even  more  pronounced  for  k  =  5. 
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TABLE  1 .  Values  of  +aE(v^  {  ))  with  k  =  3,  w  =  e^,  t  =  4.578 

and 


U1  =  (2,1,2)  U2  =  (1.5, 1,2) 


a 

+ 

p 

m 

1  H,) 

A 

n12 

v2  +  aE(v2  |  H-j ) 

A 

n12 

-1 

(.6124,-6124 

,0) 

.866 

.1448 

(.3536, -.3536,0) 

.500 

.0880 

0 

(2.449,1.225 

,2.449) 

1 

.000 

.1338 

(2.121,1.414,2.828) 

1.000 

.0481 

1 

(4.287,3.062 

,4.899) 

1 

.323 

.0770 

(3.889,3.182,5.657) 

1.803 

.0272 

2 

cn 

no 

-P* 

CO 

to 

«o 

,7.348) 

1 

.732 

.0466 

(5.657,4.950,8.485) 

2.646 

.0227 

5 

(11.64,10.41 

,14.70) 

3 

.123 

.0270 

(10.96,10.25,16.97) 

5.220 

.0192 

10 

(20.82,19.60 

,26.94) 

5 

.568 

.0213 

(19.80,19.09,31.11) 

9.539 

.0178 

TABLE  2.  Values  of  Ui2(v  +  aE(vi  H-j))  with  k  =  5,  w  =  e5,  t  =  7.665 
and  v  =  n/A(n) . 

U  =  (3, 1,3, 4, 5) 


a 

A 

tt12 

-1 

(.3371,-3371,0,0,0) 

.48 

.0699 

0 

(1.011, .3371, 1.011, 1.348, 1.686) 

1.00 

.0334 

1 

(1.686,1.011,2.023,2.700,3.371) 

1.82 

.0204 

2 

(2.360,1.686,3.034,4.045,5.057) 

2.68 

.0151 

5 

(4.384,3.708,6.068,8.090,10.11) 

5.30 

.0108 

10 

(7.753,7.079,11.12,14.83,18.54) 

9.68 

.0101 
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4.  CONTRAST  TESTS.  Suppose  one  is  to  test  HQ  versus  H^Hq  with  a 

L  _ 

contrast  test  which  rejects  for  large  values  of  T  =  L.  w.c.X.  with 

C  1  ”  t  ITT 

2 

wi  =  n^/a^  and  C  t  0.  Assuming  the  weights,  w..,  are  equal,  Abelson  and 
Tukey  (1963)  found  that  the  optimal  contrast  coefficients  are 
c}°)  «  ((i-l)(k-i+l))^^-  (i(k-i))^,  1  s  i  ^  k.  Schaafsma  and  Smid 

(1966)  generalized  their  work  to  the  case  of  unequal  weights  and  obtained 

(4.1) 


with  s.  =  Ej  =  -|  w.  and  Sq  =  0. 

We  note  that  w.cj^  =  0  and  so  the  distribution  of  T  /q^  is  the 

c' 

same  for  all  n  €  Hq. 

One  could  also  consider  testing  versus  by  rejecting  for 
large  values  of  such  a  statistic.  Of  course  the  contrast  coefficients  for 
testing  versus  H?  would  be  different  than  those  chosen  for  testing 

Hq  versus  H^-Hq.  The  power  function  for  the  test  (whether  testing  Hq 
versus  H^-Hq  or  H-j  versus  H£)  is  given  by 

(4.2)  nc(n)  -  l-4((t-(c,u)w)/|lc|Iw). 


Since  the  distribution  of  Tc  may  not  be  the  same  for  all  u  €  H^ ,  the 
level  of  significance  is  sup^H  ttc(}j).  I'f  there  is  a  u  €  H^  with 

(c,n)^  >  0,  then  using  the  fact  that  H^  is  a  cone,  we  see  that  this 
supremum  is  1.  Thus,  we  restrict  attention  to  c  with  (iu,c)^  £  0  for 
all  u  €  H.j ,  or  equivalently  c  f  Hp  For  such  c,  the  level  of  sig¬ 
nificance  is  sup^  [1  -  *((t-(c,n)w)/l|c||w)}  =  1  -  ^(t/[[c|lw).  Thus,  if  zp 

satisfies  $(z  )  =  1-p,  then  t  =  z  JjcjL  gives  a  test  of  size  p. 

P  P 


We  now  consider  different  optimality  criteria  and  the  corresponding  c. 


Fix  n  £  Hi  and  consider  the  contrast  test  which  maximizes  the  power  at  n, 
that  is,,  c  maximizes  (c,n)w/||cj[w  over  all  c  €  H*W -  { 0} .  Since  n  £  H-j , 


there  eixsts 


Consider 


i  f  j,j+l  and  c^.  =  -cj+]  =  then  p(c,u)  >  0.  If  we  agree  that 

*W 

(0,n)w/||0{[w  =  0,  then  the  maximization  problem  is  unchanged  if  H-j  -  {0} 

.  i-  U  /v 

is  replaced  by  .  Since  [j,  is  fixed  and  c  =  0  for  c  €  ,  the 


above  is  equivalent  to 


(4.3)  maximize  p(c,u)  =  w-j  (ci -o')  (ni -pT)/ ( li c-c’|[ w[[u-JjT[| w)  with  c  €  H*W 

(set  p(0,n)  =  0).  Clearly,  H*W  ©  HQ,  which  is  characterized  in  (3.3), 
is  a  closed  convex  cone  containing  the  constant  functions.  Furthermore,  c 
solves  maximize  p(c,pi)  over  c  €  H*W  ©  HQ  if  and  only  if  c-cf  solves 

(4.3) .  Applying  (ii)  of  Corollary  E,  p.  320  of  Barlow  et  al.  (1972), 

Ew(u  !  H*W  ©  Hq)  maximizes  p(c,pi)  for  c  €  H*W  ©  HQ.  Using  (2.1)  it  is 
easily  shown  that  E^(n  j  H*^  ©  Hq)  =  E^(pi  |  H*^)+iI.  Since 

2k=1  wi(Ew(u  |  H*W)i  +  PT)/E^=1  Wi  =  u,  Ew(ia  [  H*W)  solves  (4.3).  The  power 
function  of  the  resulting  test  is  1  -  *(zp-(Ew(uI  H*W)w,u)/||Ew(u  |  H*W)||W), 
which  by  (2.2)  can  be  written  as  1-  §(zp-J[Ew(u  [  H*W)[|W).  We  have  proved: 

Theorem  4.1.  Let  (j  £  H^.  The  contrast  test  with  maximum  power  at 
u  is  determined  by  c  =  E^(n  |  ).  The  power  function  is 

nc(u)  =  l-4(zp-||Ew(u|  H*“)||w). 

Since  the  optimum  c  depends  on  the  unknown  u,  one  could  estimate 
c  using  EW()T  j  H*W)  =  X"-EW(X  [  H1 ) .  However,  Ek=1  w..  (X.  -  EW(X  [  H-j  )i  )X. 

=  "X-  EW(X  (  =  T-j ^  (cf.  (2.1)).  Thus,  T]2  is  an  adaptive  contrast 


Next,  we  consider  the  criterion  used  by  Abelson  and  Tukey  (1963),  that 
is,  we  fix  6  >  0  and  seek  contrast  coefficients  which  maximize  the  mini¬ 
mum  power  over  points  at  distance  6  from  the  null  hypothesis,  .  So, 
we  wish  to  solve 

SUPc€Hf-tO)  0  -  ^p-(c,,,y  ||c||w)}. 

However,  we  will  show  that  for  c  €  H-j  -  {0]  and  6  >  0, 

inf  i1  c  i  !u  mi  s  tt  (u)  =  0  so  that  this  criterion  is  not  useful, 

u :  liU- E  w  (n  I H 1 )  I|  w=6  cvh/ 

Lemma  4.2.  If  c  €  H*W-  {0},  6  >  0  and  k  >  2,  then  there  exists 
3  U  £  H-|  with  [[pi  -  E^(pi  1  )||yj  “  6  ^nd  (c,Ey(pi  [  Hi )  )^  <  0. 

Proof.  Let  v]  =  (w^1 ,-w"1 ,0, • * • ,0) ,  v2  =  (0,w'] ,-w"1 ,0, • • • ,0) , 
■">vk-i  =  (0*  ** ,-w^).  It  is  easy  to  show  that 
H*W  =  {aivi  +  •••  +  ak_ivk_i  :  ai  21  °i-  Furthermore,  (v.  ,Ew(n  {  H-,  ))w  s  0 

for  each  i  and  pi.  Let  c  =  a^  +  ’ '  ’  +  ak  i  v(<  l  *  ^  aj  >  0> 

1  £  j  <  k-1,  then  let  pi  be  (1,2,  *•  *,k)  with  the  j+1^  and  j+2~ 
coordinates  interchanged.  So,  Ew(ja  (  =  i  for  i  f  j+l,j+2  and 

Ew(ulHi)i  -  ((j+2)wJ+i+(j+1)wj+2)/(wj+i  +  wj+2)  for  i  =  j+l,j+2.  Thus, 

a j ^ v j l  ^ =  aj ( j-E^(n  [  Hi ) j+i )  <  0  and  (c,E^([i|Hi))^<0.  If 

ai  =  •••  =  ak_2  =  0,  and  ak_i  >  0  (recall,  c  f  0),  then  let 
u  =  (2,1,  —  ,k).  It  is  easy  to  show  that  ak_-j(vk_i,Ew(ial  H1))w  <  0. 

Thus,  in  either  case,  one  can  find  pi  £  Hi  with  (c,E^(pt[  H^y  <  0. 
Multiplying  by  the  appropriate  positive  constant,  we  obtain  pi  £  H.^  with 
(c,Ew(pt|  Hi))w  <  0  and  ||ja-Ew(n  {  )|JW  =  6.  The  proof  is  completed. 


For  a  £  -1,  set  pa  =  p  +  aE^(p[  H-j )  and  note  that  by  the  remark  pre¬ 
ceding  Theorem  3.1,  l^a"EW^a  ^  H1^W  =  ll^_Ew^  l  Hi  Miw  =  6  >  °-  Thus» 
u  €  H-, ,  the  distance  from  p  to  H.  is  6  and 

a  I  a  I 


=  (c"^w  +  ,1m^a(c-Ew(|Jl  Hl>>w  =  - 


Therefore,  for  each  c  €  H*W- £ 0}  and  6  >  0, 


V:  lil^-Ew(ul Hn ) Ilw-6  c 


—  c  tt  (y.)  =  0. 


We  must  consider  other  criteria. 


Following  Schaafsma  and  Smid  (1966),  we  consider  the  contrast  that  mini 


mizes  the  maximum  "shortcoming"  among  all  contrast  tests.  Recall  that  for  a 


given  p  (6  H-| ,  the  contrast  test  with  maximum  power  at  p  is  obtained  by 
taking  c  =  p-Ew(p  |  )  and  has  power  1  -  §(z  -||p-Ew(p  |  H-j ) Uw)  •  So,  for 
any  contrast  test  its  shortcoming  at  p  is 


(4.4)  $(z  -(c,p)w/|Ic!!w)  -  Hz -||p-Ew(p  I  H1)||w). 


If  there  is  no  constraint  on  p  other  than  p  ,  we  see  from  the  pre¬ 


ceding  analysis  that  the  supremum  is  at  least  as  large  as  1  -  $(zp-&)  for 


each  6  >  0,  and  so  the  maximum  shortcoming  over  all  p  £  is  1.  Even 
if  p  is  constrained  so  that  [jp-E^(p  [  H^)(|^  =  6  >  0,  the  maximum  short¬ 


coming  is  1  -  §(zp-6)  which  does  not  depend  on  c.  Neither  of  these  cri¬ 


teria  are  useful 


The  vector  of  means  p+aE^(p  [  )  remains  at  a  fixed  distance  from 

,  but  it  is  moving  away  from  as  a  increases.  So,  we  consider  the 


contrast  test  which  maximizes  the  minimum  power  over  all  p  £  with 


A(p)  =  |[p-pjjw  =  6  >  0.  Let  ai  =  (w^+wT^  )^2  for  i  =  1 , 2 ,  *  *  * , k-1 


.  _  .  %  -  •  •  •  - 
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let  d^  =  0,  =  Lj~]  for  i  =  2,---,k,  let  d^  =  d'-cf'  and  let 

c^.-d,. 

Theorem  4.3.  Let  6  >  0.  The  contrast  test  which  has  coefficients 

c^  and  rejects  for  large  values  of  T  /■,*  maximizes  the  minimum  power 

c '  ' 

over  all  n  £  H-|  with  A(pi)  =  6.  Furthermore,  such  contrast  coefficients 
are  unique  up  to  a  positive  multiplier. 

Before  the  proof  of  the  theorem  is  given,  we  establish 

Lemma  4.4.  If  pi,v  €  H^,  then  (pi-u,v-v)w  ^  0. 

~  *W  ~ 

Proof.  n-|a  6  and  v-v  €  and  the  conclusion  is  immediate. 

Proof  of  Theorem  4.3.  We  wish  to  find  c  which  yields  sup 

-  -  c€H*W-£0] 

inVH1,A(n)=6^"i^zp_(c’w^W/HcIlw^’  or  since  c  =  0,  equivalently, 

(4.5)  sup  w  inf  p(c,u). 
c€H*W-£0j  ^H1 

If  -c  f.  ,  then  consider  pi  =  -c.  Since  p(c,-c)  =  -1,  we  may  omit 
such  c  from  the  supremum.  Because  H*W  0  (-H^)  =  (-H^)  fl  £u:|J=0}, 

(4.5)  is  equivalent  to 

(4.6)  suP-c€H1-£0},c=0  infu^H1  p(c’^  =  -1  nfd€H1  -£0} ,3=0  SUpia)6H1 

We  will  solve  for  d  and  remember  that  c  =  -d.  Because  of  the  continuity 
of  p(d,*)»  the  supremum  in  the  r.h.s.  could  also  be  taken  over  u  (6  H^UH^. 
However,  if  n  €  ,  then  applying  Lemma  4.4,  p(d,u)  s  0  and  so  that 

supremum  could  be  restricted  to  pi  £  U  (Hq  ©  H*W)  (for  p(d,u)  s  0 
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*W 

for  p,  6  ).  So,  we  seek  d  which  solves 

(4.7)  inf .  „  fn-j  ^_n  sup  n  u  p(d,n). 

d€Hr10j,d-0  ^h°U(H0G«*W} 

Furthermore,  if  u  €  H*W  0  HQ,  then  Ew(u  {  H^)  f.  HQ  ar.d  so  A(Ew(n  (  H-j ) ) 

k  k 

>  0.  Applying  (2.3)  and  the  fact  that  E.  ,  w . E, . (n  [  H,).  =  £._>  w.pi.,  we 
see  that  0  <  A(Ew(u  |  H-j))  =  ||Ew(n  |  H^-juUy  =  l|Ew(u“Ct  I  H-j )  11  w  ^  11^'Slw  = 

For  fixed  d  €  H^-{0}  with  cT  =  0,  p(d,u)  =  ( I|d|| w •  I;u-ul!w) ” 1  (d,n)w 
£  (||djjw  lin-il^)"1  (d,Ew(n  [  H^))w>  which  is  nonnegative  by  Lemma  4.4.  So 
p(d,u)  £  p(d,Ew(n(  H-j ))  for  n  j£  hQ.  Therefore,  d  solves 

infd€HrC0}>0  *“Plia^(H^)P<d'EW<,‘,Hl» 

=  Tnfd€H1  -{03 ,H=0  SU^€dH1-HQ  p(d’^’ 

The  boundary  of  H-j  is  the  union  of  A-j  =  [x  €  R  :  x-j  =  x2 ^  x3 s  •  •  •  s;  xkJ, 

k  r  k 

A2  =  [x  €  R  :  x]sx2=  x3^  •••<:xk3,'--,Ak_1  =  IX  €R  :  X]  £  .  •  •  s  x|<_1  =  xkJ . 
Because  of  the  convention  p(d,0)  =  0,  we  seek  d  that  solves 

(4.8)  1<*()€hi.£oJ,3=o  maxUisk-l  o(d’u)' 

k  • 

Each  A^  is  a  closed,  convex  cone  in  R  containing  the  constant  func¬ 
tions  and  p(d,j_i)  ^  0  for  any  n  €  A. .  So,  by  Corollary  E,  p.  320  of 
Barlow  et  al.  (1972),  max^^  p(d,u)  =  p(d,E^(d  |  A^ ) ) .  It  is  easy  to 

show  that  d*  =  (dt\ •  •  •  ,d£) ,  with  d1*  =  d.  for  j  t  i,i+l  and 
IK  J  J 

dj  =  (widi+wi+idi+1 )/(wi+wi+i )  for  j  =  i  ,i+l ,  is  the  point  in  A. 
closest  to  d  €  ,  i .e. ,  d*  =  E^(d  [  A.j ) .  Also 


p(d,Ew(d  |  Ai))  =  }|Ew(d  |  Ai)l|w/||d|lw.  So  d  solves  (4.8)  if  and  only  if 
d/i|d|iw  solves 

""oeH,  ,Ud|iM-i .3=0  ",axiSi£k-iUEu(d  I  Ai >lltr 

However,  ||d||w  -  p„(d|  At)Bw  =  w.w1+] (d1+1-d,)2/l*(1+wltl)  -  @j+,-d,)2/a?. 
So,  we  wish  to  solve 

(4-9>  maxdeH,  ,|idSw=i  .3=0  mini*isk-i(di+rdi>2/ai- 


Let  d]  be  defined  as  in  the  paragraph  before  the  statement  of  the  theorem 
and  let  dg  =  ad-].  Note  that  dfl  €  ,  cfa  =  0  and  |jda||w  =  a  •  ]Jd]J|w  >  0 
for  all  a  >  0.  So,  choose  a  so  that  J|d  ||^  =  1. 

We  now  show  that  the  d„  chosen  above  is  the  unique  solution  to  (4.9), 

a 

which  implies  that  -da  is  the  unique,  up  to  a  positive  multiplier,  set  of 
contrast  coefficients  which  is  being  sought.  Note  that  if  dg  =  (dai »da2 ’ 

•  *  * ,dflk) *  then  (dfll.+1  -  da_.)2/a2  =  a2  for  i  =  1,2,-*- ,k-l .  Suppose 
z  €  H]  with  z  =  0,  |[z|lw  =  1  and  minl£.i^k_-| (z.+1-zi )  /a.  *  a  .  Then, 

(2i+rzi>  *  <dai+rd«i>  or  zi+rdai+i  *  zrd«,  for  i  ’ 

Hence,  z-d,  €  H-.  and  applying  Lemma  4.4, 

a  I 


1  =  Kw  ’  llda»U  + 


+  2(da,z-da)w  *  1+ilz-dal(w- 


So,  i|z-daiiy  =  0  or  z  =  da. 


The  proof  is  completed. 


We  conclude  this  section  with  some  remarks  concerning  the  power  func¬ 
tions  of  such  contrast  tests,  that  is,  tests  which  reject  for  large  values 
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i/ 

Theorem  4.5.  Let  pi,v  €  R  .  If  (c,v)u  =  0,  then  it  (pi  +  av)  =  tt  (n) 

w  c  c 

for  all  a  If  (c,v)w  >  0  ((c,v)w<0),  then  TTc(u  +  av)  is 

increasing  (decreasing)  in  a  with  lim  tt  (pi  +  av)  =  1  (0)  and 

TTc^  +  av)  =  0  O)- 

Proof.  The  result  follows  immediately  from  (4.2)  since  (c,u  +  av)y 
=  (c,u)w+a(c,v)w. 

In  the  next  result  the  regions  of  consistency  are  determined  for  such 
contrast  tests. 

Theorem  4.6.  Let  n,y  €  with  y.  >  0  for  i  =  l,2,***,k.  Let 

wn  =  r)  y  and  fix  the  level  of  the  contrast  test  at  p  €  (0,1)  for  all  n. 

If  (c,n)  >  0  ((c,n)  <  0),  then  n  (pi)  1  (0)  as  n ®.  If 

Y  Y  c 

(c,)j)  =  0,  then  tt  (n)  =  p  for  all  n. 

Y  c 

1/2 

Proof.  Since  (c,ja)w  /||cjL  =  n  (c,pi)  /||cj|  ,  the  desired  conclusion 

n  n  Y  Y 

follows  from  (4.2). 

It  is  of  interest  to  compare  the  regions  of  consistency  for  T^  and 
T  ^  in  testing  Hq  versus  H-j -H^ .  We  first  show  that  c^  €  H^.  Let 
x^  =  s^/s^  for  i  =  0,1,  •••,!<  and  note  that  (cf.  (4.1)) 

S(0)  «  (g(xi_i)-g(xi))/(xi-x._1)  where  g(x)  =  (x(l-x))1^2.  The  desired 
conclusion  follows  from  the  strict  concavity  of  g  on  [0,1].  Applying 
Theorem  4.6  the  contrast  test  is  consistent  for  pi  €  A  ic^) 

=  Cm  *  (</°\u)w  >  0],  and  from  Corollary  3.2,  thi^  is  true  for  T^  for 
U  £  H*W  ©  Hq.  If  pi  =  pt'  +  u"  with  pi'  €  H*W  and  pi"  €  HQ,  then 

*  (</°\pi#  )w  £  0.  So,  A+(c^)  c  (H*W©Hq)c.  Drawing 


H1’  »  c  and  A  (c  )  for  k  =  3  in  the  plane  jJ  =  0  gives  one 

an  idea  of  the  size  of  (H*W©HQ)C  -  A+(c^ ) . 

In  testing  H-j  versus  H2,  "*"(!)  consistent  for  n  6  A+(c^) 

c 

and  T]2  is  consistent  for  w  £  H] .  Since  - c ^  €  H]  and  c^  =  0, 
£1}€H*W.  If  u€Hr  then  (c(1)4i)w<:0  and  n  e  (A+(c(1  * ) )c.  Henc 
A  (c^)  c  hJ.  Again  one  can  obtain  an  idea  of  the  size  of  H^-A+(c^) 
by  drawing  the  figure  for  k  =  3. 
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5.  COMMENTS.  We  begin  this  section  with  a  few  remarks  about  the  situa¬ 
tion  in  which  the  variances  are  unknown.  Suppose  X.  .  are  independent 

*  J 

2 

r)(|j.j,CT  )  variables  for  j  =  l,2,***,n.  and  i  =  l,2,*--,k  with 

p 

H  =  (u-j ,U2> • • • »uk)  and  a  unknown.  For  the  contrast  tests,  let  w^  =  n.. 

2 

and  assume  that  c  does  not  depend  on  a  .  Following  the  optimal  procedure 
for  k  =  2,  define  c2  =  £k  (X.  .-X. )2/(N-k)  with  N  =  £k  n . .  In 

^  ^  v  "I  Til 

b 

testing  HQ  versus  H.|-H  ,  we  assume  Z-_i  W..C..  =  0,  and  so  if  one 

lr  _ 

rejects  for  large  =  Z^_-j  w^c^X./5,  then  the  1 00p%  critical  value  is 

t..  .  jjcIL  where  F(t..  .  )  =  1-p  with  F  the  distribution  function  for 

N-k,pM  "W  N-k,p 

a  Students  t  variable  with  N-k  degrees  of  freedom.  In  testing 

*W 

versus  H2,  we  assume  c  €  H-|  ,  and  so  Hq  is  least  favorable  within  . 

Hence,  the  100p%  critical  value  is  also  t^_^  ||c|ly.  Let  f(y)  be  the 

density  of  Yn  =  ct/ct.  Conditioning  on  YN>  which  is  independent  of 
X  -  (Xyl2,^  •  •  ,x'k) ,  we  see  that  the  modified  contrast  test  has  power  func¬ 
tion 


(5.1)  nj.(u)  =  1  -  f  i^tN-k,p"  (c»u)w/^c[[w))f(y)dy. 


Hence,  Theorem  4.5,  which  gives  the  radial  monotonicity  and  radial  limits 

of  ttc,  is  also  valid  for  t/.  Furthermore,  if  n^  =  ny^  with  y-j  >  0 

for  i  =  1 ,2, •  •  •  ,k,  then  as  n  ->  *,  £k=1  rucX/Hcl^  ■  Jr\  z{=1  YiciX./i!c[!Y 

->  ±°°  depending  on  whether  (c,n)^  >  0  or  (c,n)^  <  0.  Also 

(a/cr)tN_k  p  ->  zp.  So  as  n  ®,  n^(ja)  -s>  1  (0)  as  n  -$»  ®  provided 

(c,n)  >  0  ( (c,ui)  <  0),  and  tt*  (u)  p  if  (c,n)  =  0.  The  radial 

Y  Y  c  Y 

behavior  and  the  regions  of  consistency  for  these  modified  contrast  tests 
are  like  those  for  the  contrast  tests. 

The  LRT  for  Hq  versus  H-j-Hq  rejects  for  large  values  of 


*  >  v 


••  -  '>.*•.•.•  .  ■  ■  .  V  ■  V  > 


S01  ~  ^i=l  ^ j = 1  ( X i j "^ ) 2  =  a?T01/ST 

where  Sj  is  the  total  sum  of  squares  (cf.  Barlow  et  al.  (1972),  p.  121)). 
Equivalently,  one  could  reject  for  large  values  of  LQ1  =  (N-k)S01/(l-S01 ) . 
But  ST  =  Ej=1  Zj^1(Xij-u.)2  +  E^=1  ni(ni-ji)2+  2Z^_i  n1(pii-)3)(Xi-u1).  Apply¬ 
ing  (2.1)  and  the  fact  that  Z^-|  w.Ew(x[  c)..  =  zJjL-j  w.x..,  we  see  that  the 
last  sum  is  zero.  Hence,  LQ1  =  (N-k)Z^=1  n.  (jl. -ja)2/!^  (X.  .-jl.  )2.  We 

now  wish  to  determine  the  region  of  consistency  for  .  Suppose  that 
ni  =  ny^.  with  >  0.  Noting  that  Z*=1  Ejl-j (X^-j^ r/(N-k) 

=  Z^=1  Zjl1(X..-Xi)2/(N-k)  +  Z^1  n.  (X\-jI.  )2/(N-k) ,  we  see  that  this  expres- 

2  2  —»k 

sion  converges  almost  surely  to  a  +|jja-E  (u|  H-j)]^/- (=1  y^ .  Under  HQ 

(we  may  assume  without  loss  of  generality  that  |j  =  0),  Z^  n^ (JT^ -|a)2 

=  Y1(EY(vlHi)j-ZJk,1  Yj)2  where  V,  X,,- •  •  ,Vk  are  inde- 

pendent  variables  with  Y..  ~‘n(0,y;:1).  Hence,  the  100p%  critical  value 

for  Lq.|  converges  to  the  1 00p%  critical  value  for  TQ1  with  weights  y.. . 

*W 

Now  suppose  n  €  ©  Hq.  Examining  the  proof  given  in  Barlow  et  al . 

(1972)  for  their  Theorem  3.4,  we  see  that  Z^=1  n.(jL-ji))2  -a‘--'  >  ®  and, 

1111 


®  and. 


as  we  have  seen,  Z-=1  (X.^-ia.  )2/(N-k)  — -S-‘->  o2+  |[u-E y(u  I  H]  )||^/E^=1  Yj 

Hence,  is  consistent  for  such  pi.  Furthermore,  the  argument  given  in 

the  second  part  of  their  proof  shows  that  Z^_-|  n_.(|j. -£)  -  ■ ---  >  o  pro¬ 
vided  n  €  (H*W  ©  Hq)^.  ^his  same  behavior  as  was  observed 


for  T, 


It  is  interesting  to  note  that 


L01  '  Q+  T12 

where  Q  is  independent  of  and  Recall  from  the  introduction 


is  isotonic  and  T 


that  T 


01 


lows  that  for  fixed  q,  P 


12  is  anti  tonic  with  respect  to  <.  It  fol- 


ia[q+T 


12 


t  is  isotonic  with  respect  to  < 


and  by  conditioning  that  P^[Lqi  >  t]  is  isotonic  with  respect  to  <. 


Thus,  if  n  €  Hi  then  P^CLqi  >  t]  is  nondecreasing  for  6  €  (-<*>,“) 
and  )^Lqi  >  t]  is  nondecreasing  for  6  €  (-00,00)  for  any 


U 


€  R^.  What  about  directions  like  n+6Ew(n  |  H*)  or  for  p.  €  (-H*)  U  Hn? 


'1 


‘1 


It  is  easy  to  see  that  for  fixed  q  >  0,  {  x  €  R  : 


k  ‘‘EW^C01^W 


9  +  ?)ll 


£  t 


12y|lW 


} 


is  not  convex  (take  q  =  t  =  1  and  k  =  3)  so  that  the  techniques  of  Sec 
tion  3  will  not  apply. 

The  LRT  for  Hi  versus  H 2  rejects  for  large  values  of 

si2 =  "i<v^)2/[(N-k>s2+4=i 


(cf.  Robertson  and  Wegman  (1978)),  or  equivalently  for  large  values  of 

I_i2  =  (N-kJS^/d-S^)  =  lki  =  ]  n.(X.-p.)2/a2  =  o2T ^/a2. 

If  we  again  denote  the  density  of  =  S/cr,  by  f(y),  then,  for  a  fixed 
critical  value,  t  >  0,  the  power  function  of  Ug  is  given  by 

TT] 2 (p )  =  >  ^ 

Hence,  Theorem  3.7,  Corollary  3.8  and  Corollary  2.11  are  also  valid  for 
rr'^.  **  not  dl^1cu^  t0  show  that  Si g  is  consistent  for  all  pi  €  Hi 

as  was  the  case  for  T^. 

Because  of  the  similarities  we  have  observed  between  the  case  of 
variances  known  and  the  case  of  variances  unknown,  one  might  conjecture 
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that  the  monotonicity  properties  of  are  like  those  of  .  However, 

we  have  seen  that  some  of  the  results  do  not  follow  from  simple  condition¬ 
ing  arguments  as  in  the  case  of  and  L^.  It  would  be  of  interest  to 

know  what  techniques  could  be  applied  in  the  study  of  the  monotonicity  of 
the  power  function  of  Lq^ . 

In  deriving  the  optimal  contrast  test  for  versus  H^,  the  vector 
d-j  was  obtained.  This  vector,  which  in  the  case  of  equal  weights  has  uni¬ 
form  increments  (i.e.,  d-j  i>+^-d^  .  is  constant),  is  in  the  center  of  the 
cone  .  In  fact,  d-j  makes  equal  angles  with  the  faces  of  .  On  the 
other  hand,  the  optimal  conrast  test  of  Hq  versus  is  based  on 

c(Qy  which  is  another  center  of  .  The  vector  c^j  makes  equal 
angles  with  the  edges  of  H-| .  Bartholomew  (1961)  conjectured  that,  for 
a  fixed  value  of  A,  the  power  of  Tq-j  is  largest  at  d^ .  It  is  of 
interest  to  compare  the  power  of  T^  at  both  of  the  "centers"  mentioned 
above.  Fixing  their  lengths  to  be  1  and  k  =  5,  d^  =  (-.6325, -.3162,0, 
.3162, .6324)  and  c^  =  (-.6899, -.1551  ,0,.  1551  ,.6899).  For  w  =  e5> 
these  powers  were  estimated  by  a  Monte  Carlo  experiment  with  9,999  repli¬ 
cations.  The  estimates  are  nQ-j(d^)  =  .2374  and  ^oi ^c(0) ^  = 
which  tends  to  confirm  Bartholomew's  conjecture.  What  analytic  tools  are 
needed  to  establish  this  conjecture? 


APPENDIX.  The  appendix  contains  the  proofs  that  were  omitted  in  Sections 
2  and  3. 


Proof  of  (2.4).  It  follows  immediately  from  the  definition  of  a  dual 
cone  that  C  c  (C  )  .  The  other  containment  depends  on  the  fact  that  C 
is  closed.  Suppose  x  €  (C*)*  and  x  £  C.  Since  C  is  closed, 

(I x-E (x  [  C)|[  >  0.  But,  x  €  (C*)*  and  x-E(x  [  C)  €  C*  imply  that 

0  *  (x,x-E(x  |  0)  =  jjx-E(x  [  C)H2+  (E(x  |  C),x-E(x  [  C))  =  ||x-E(x  I  C)|j2.  Thi 

contradiction  shows  that  (C*)*  c  C. 

Proof  of  Lemma  2.2.  We  first  note  that  if  E(x-(Jq  [  S  )  =  hpi,  then 
E(x-|aQ[  C^)  =  b'u  where  b'  =  bvO,  and  E(x-hq  [  S  )  =  E(x  |  S^).  Hence, 

E(x-|_igl  C  )  =  E(x  |  C  )  and  so  we  establish  (2.5)  with  =  0. 

We  consider  the  two  cases  n  €  C  and  -pi  €  C*  separately.  Suppose 
u  €  C  and  0  fi  b  s  2.  Using  (2.3)  followed  by  Lemma  2.1  and  (2.3)  again, 
we  see  that 


|jE(x-bE(x  (  y  [  C)lr  =  ||x-bE(x[  C^)||  -  |[E(x-bE(x  l  C^)  I  C*)|f 

*  i|x-bE(x{  Cu)||2-  ||E(x[  C*)||2 
=  IJxjj2  +  b(b-2)[|E(x  1  C^)|)2  -  U  E  ( x  [  C*)||2 


xf-  j|E(x[  C*)||2  =  |!E(x|  C)||2. 


If  u  €  C*,  then  -bE(x  [  C  )  €  C*  for  all  b  £  0.  Thus,  by  Lemma  2.1 


and  (2.4), 


HE(x-bE(x[  y  (  C)l|  *  HE(x  i  C)U 


for  all  b  £  0. 
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Proof  of  Lemma  2.3.  Because  of  (2.3),  ||E(x+y[  C)[[  can  be  written  as 

(A.  1 )  ||E(E(x  |  C)+E(x  {  C*)+E(y  [  C)  +  E(y  (  C*)  [  C)|| 

=  ii E ( E (x  [  C)+E(y[  C)  +  z  f  C)|i 

where  z  =  E(x  [  C*)  +  E(y  [  C*)  6  C*.  Applying  (2.4)  and  Lemma  2.1,  (A.l )  is 
bounded  above  by  J|E(E(x  [  C)  +  E(y  [  C)  {  C)|j  =  |[E(x  [  C)  +  E(y  |  C) || .  The  second 
inequality  in  Lemma  2.3  follows  from  the  triangular  inequality  for  norms. 

Proof  of  Lemma  2.4.  The  first  conclusion  of  part  (a)  follows  from  the 
third  condition  in  (2.1)  and  the  facts  that  -v  €  S  whenever  v  €  S  and 
ScC.  For  the  second  conclusion  in  part  (a),  we  check  the  three  conditions 
in  (2.1).  Clearly,  E(x  (  C)-v  €  C  and  (x-v-(E(x  {  C)-v) ,E(x  {  C)-v) 

=  (x-E(x  {  C) ,E(x  (  C))-  (x-E(x  (  C),v),  where  the  first  term  on  the  r.h.s 
is  zero  by  (2.1)  and  the  second  term  is  zero  because  of  the  first  part  of 
(a). 

For  part  (b),  we  assume  S  is  a  closed  subspace  contained  in  C  and 
show  that  E(x  (  S)  satisfies  the  three  conditions  that  characterize  the 

projection  of  E(x  {  C)  onto  S.  Of  course,  E(x  [  S)  €  S, 

(E(x  |  C)-E(x  |  S) ,E(x  |  S))  =  (x-E(x  (  S),E(x(  S))-(x-E(x[  C),E(x|  S))  =  0 
(recall ,  E(x  (  S) ,-E(x  [  S)  €  S  c  C) ,  and  for  u  €  S,  (E(x  [  C)-E(x  [  S) ,u) 

=  (x-E(x  I  S),u)-(x-E(x[  C),u)  =  0  (again,  u,-u  €  S  c  C) . 

We  prove  part  (d)  before  (c).  So,  we  assume  that  CcS  and  again 
verify  the  conditions  in  (2.1).  By  definition,  E(x  [  C)  €  C,  and  because 
x-E(x j  S)  =  E(x  (  S1) ,  (E(x  I  S)-E(x  (  C),E(x  {  C))  =  (x-E(x  [  C),E(x  {  C)) 

-  (x-E(x  [  S),E(x  [  O)  =  -(E(x  J  SJ’),E(x  |  O)  =  0  since  CcS.  For  y  €  C, 

( E ( x  |  S)-E(xf  C),y)  =  (x-E (x  (  C),y)-  (E(x |  Sx),y)  =  (x-E(x  [  C),y)  *  0. 


For  part  (c),  assume  that  S  is  a  closed  subspace  contained  in  C. 

By  part  (a) ,  E(x  (  C)-E(x  [  S)  =  E(x-E(x  [  S)  {  C)  =  E(E(x  {  S1)  [  C) .  By  part 
(a),  for  v  €  S,  (E(E(x  |  SA)  |  C),v)  =  (E(x  |  Sx),v)  =  0  and  so , 

E(E(x  |  SA)  (  C)  €  C  H  Sx.  Hence,  E(E(x  [  Sx)  [  C)  =  E(E(x  I  Sx)  |  C  n  Sx)  and 
the  latter  is  E (x  {  C  fl  Sx)  by  part  (d). 

Proof  of  Lemma  2.9.  For  the  proof  of  the  first  conclusion,  we  note 

that  x  e(FC)*1  (x,y)  £  0  for  all  y  €  FC  <=>  (x,Fz)  £  0  for  all 

-1  +W 

z  €  C  <=MF  x,z)w  sO  for  all  z  6  C  <=>  x  €  FC  . 

For  the  second  conclusion,  we  show  that  FE^(x  (  C)  satisfies  the  three 

conditions  that  characterize  E(Fx|  FC)  (cf.  (2.1)).  Let  y  =  Ew(x  {  C). 

Of  course,  Fy  6  FC,  (Fx-Fy,Fy)  =  (x-y,y)w  =  0.  and  for  z  €  FC, 

(Fx-Fy,z)  =  (x-y,F-1z)u  s  0  since  F-1z  €  C. 
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